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ABSTRACT 


Several  programs  are  *n  existence  which  perform  quantal  analysis, 
but  such  programs  were  designed  for  batch  processing.  The  following 
report  is  a  package  of  computer  programs  for  the  IBM  2250,  with  documen¬ 
tation,  which  enables  the  layman  to  interact  with  the  system.  The 
program,  entitled  "QUANTAL"  is  written  in  FORTRAN  IV  with  extensive  use 
of  the  GRAF  (Graphics  Addition  to  FORTRAN)  routines.  QUANTAL  is  written 
in  subroutine  form  within  an  overlay  structure  which  retains  the  entire 
program  size  to  within  64K  bytes.  Communication  between  subroutines  is 
attained  by  extensive  use  of  COMMON  variables  and  direct  disk  storage. 

The  method  of  analysis  follows  the  usual  techniques  of  maximum- likelihood 
iteration,  with  an  approximate  Newton  technique. 
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CHAPTER  I 
INTRODUCTION 

i  great  variety  of  experiments,  especially  those  in  biological 
assay,  establish  a  relationship  between  dosage  (a  concomitant  variable) 
and  a  quantal  response.  The  latter  consists  of  records  of  successes  or 
failures  under  a  given  dosage.  The  quantal  response  is  often  ."death," 
but  may  be  any  other  easily  recognized  change  in  an  experimental  subjecc. 
For  example,  insecticides  may  be  assayed  by  assigning  batches  to  doses 
of  a  test  preparation  --  here  the  quantal  response  is  the  number  of 
insects  killed  in  each  batch.  In  the  physical  or  social  sciences,  t!  e 
concomitant  variable  is  usually  "time."  At  specific  intervals  a  record 
is  kept  about  the  number  of  subjects  that  completed  a  task,  or  the  number 
of  experimental  runs  that  have  proceeded  to  completion  in  a  reaction. 
These  are  the  quantal  responses. 

This  thesis  first  summarizes  a  few  of  the  well  known  basic  results 
for  the  analysis  of  quantal  response  data.  Essentially,  the  approach 
presented  by  Finney  [3],  [4]  is  followed,  including  the  notation,  some¬ 
what  peculiar  for  statisticians,  which  uses  capital  letters  for  para¬ 
meters.  Chapter  II  also  includes  a  discussion  of  the  various  algorithms 
as  they  apply  to  a  few  preferred  growth  functions  (Probit,  Logistic, 

Arc  Sine,  Log-Log,  Weibull  --  see  [1],  [3]). 

The  implementation,  including  user's  guide,  description,  and  illus¬ 
tration,  is  contained  in  Chapter  III.  This  is  an  example  of  an  interac¬ 
tive  program  designed  for  the  user  whose  expert  knowledge  in  Statistics 
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or  Computations  may  be  limited.  In  fact,  the  present  program  has  been 
used  in  the  research  work  of  entomologists  at  the  University  of  Georgia. 
This  part  of  the  work  is  directed  toward  the  implementation  of  the  Univer¬ 
sity  of  Georgia's  development  of  Statistical  Programs  for  Lay  Users,  under 
the  THEMIS  program. 

Detailed  documentation  of  the  various  computer  programs,  the  communi¬ 
cation  with  the  Graphics  Monitor  System,  and  the  overlay  structure,  are 
presented  in  the  last  chapter. 

A  statistician  desiring  to  use  the  computer  for  the  analysis  of 
experimental  data  has  a  variety  of  programs,  procedures,  and  languages 
at  his  disposal.  The  majority  of  these  programs  (BMD  [2],  GENSTAT  [9], 
OMNITAB  [6],  and  SSP  [8],  to  name  a  few)  require  some  degree  of  technical 
understanding.  Often  a  user  must  even  supply  a  complete  function  pro¬ 
gram  or  a  subroutine.  Thus,  unless  the  statistician  possesses  enough 
experience  in  programming  and  data  processing  to  prepare  his  input  data, 
he  must  seek  the  advice  of  one  with  such  experience. 

The  interactive  console  programs,  on  the  other  hand,  arc  designed 
so  that  by  merely  answering  questions  and  filling  in  blanks  a  layman, 
who  possesses  little  or  no  knowledge  in  the  fields  of  Statistics  or 
Computation,  may  perform  the  analysis  on  his  particular  data.  The 
QUANTAL  program  also  affords  the  layman  the  opportunity  to  change  his 
data  input  (by  adding  or  deleting  batches,  or  performing  some  transfor¬ 
mation  on  the  x-dosage)  and  immediately  see  the  results.  This  eliminates 
the  necessity  of  re-submission  of  the  job  deck  under  the  batch  mode  and 
thus  increases  the  overall  effectiveness  of  turnaround. 


CHAPTER  II 

COMPUTATIONAL  ALGORITHMS 


In  a  biological  assay,  a  stimulus  (for  example  a  vitamin,  a  drug, 
a  physical  force)  is  applied  to  a  subject  (for  example  an  animal,  a  plant, 
or  a  single  cell)  at  an  intensity  specified  in  units  of  concentration, 
weight,  time  or  other  appropriate  measure  and  under  environmental  condi¬ 
tions  as  carefully  controlled  as  is  practicable,  as  a  result  of  which  a 
response  is  produced  by  the  subject. 

The  derivation  of  iterative  maximum- likelihood  estimates  follows 
familiar  lines  (see,  e.g. ,  Finney  [3]).  Let  us  assume  that  there  are  k 
classes  of  subjects  and  in  is  the  size  of  class  i.  Suppose  r^  denotes 
the  number  of  "successful"  individuals  in  class  i  and  P.  denotes  the 

i 

probability  that  a  subject  in  class  i  has  made  the  transition.  The 
probability  of  exactly  r^  responding  is,  the  Binomial  Distribution  of 
probabilities 


ni!/ri!(ni  -  r^l 


r.  n.-r. 
P.1  Q.  1  1 

l 


where  Qi  *  1  -  P^.  Then  the  logarithm  of  the  likelihood  function  of 


(r^  i  r2»*  •  *rj<)  ^ 


log  L 


k  .  . 

I  Hog  (  i)  *  r. l°gP,  •  (n,  -  r,)log(l  -  P  )]  (2.1) 

1.1  Vr J 
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If  we  let  denote  the  observed  proportion  of  successes  in  class  i 
(p^  =  r^/n^) »  we  find  that 

3 log  L/3P.  =  n.(p.  -  P.)/P.(l  -  P.)  (2.2) 

As  is  customary  in  bio-assay  analysis,  we  use  capital  letters  to  denote 
expected  values  and  small  letters  to  denote  observed  values. 

The  dosage- response  analysis  proceeds  as  follows.  A  monotonic 
(growth)  function  is  defined 

Pi  =  *(Yi)  (2.3) 

which  is  frequently  a  cumulative  distribution  function.  It  is  then 
assumed  that  Y.^  follows  a  linear  regression  model 

=  a  +  0  xA,  (2.4a) 

where  x^  is  the  dosage  (or  some  convenient  transform) . 

Y.  =  *"1(P) 
l  i 

.  is  a  "percentage  point"  under  this  c.d.f.  It  is  a  linear  function  of  x^. 
Let  Zi  =  dP./dYi  (2.5) 

Let  x  denote  a  conveniently  chosen  constant  (usually  a  weighted 
mean  of  the  x^)  and  let 

U  =  a  +  fix,  so  that  Yi  =  \t  +  fi  (xi  -  x). 


(2.4b) 
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Then 

k 

a  log  L/ay  =  l  [aiog  L/ap.HdP./dY.naY./ati] 
i=l 


k 

-  I  [^(Pi  -  P.)Z.]/[P.(1  -  P.)]  (2.6a) 

and,  similarly, 

k 

3 log  L/33  =  . I  [Hi (Pi  -  P^fr.  -  x)]/[P.(l  -  P^]  (2.6b) 

Expressions  (2.6a)  and  (2.6b)  must  be  equated  to  zero  and  solved 
for  y  and  3  (involved  in  each  Y^,  hence  in  each  P^  and  Z^).  This  is 
customarily  done  by  an  approximate  Newton  iteration  where  y^ 
is  used  as  the  first  approximation  to 
If  we  define 

W.  *  Z*/[P.(1  -  P^L  (2.7) 

(2.6a)  and  (2.6b)  can  be  rewritten 

k 

3 log  L/3y  =  l  n  J  ((p  -  P.  )/Z.  ] 
i=l  1  1  1  1  x 

and 

k 

31og  L/33  ■  l  niWi[(pi  -  Pi)/Zi)(x.  -  x) 
i»l 

Now, 

,  2  k 

3  log  L/3y‘  »  l  { [3(31og  L/3y)/3P. ] [3P^/3Y. ] [3Y^/3y] 
i») 


(2.8a) 


(2.8b) 


♦  [3 (Slog  L/3y)/3Zi](3Zi/3Yi][3Yi/3y]}. 
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Now 

a(iog  L/aiO/aPi  =  -(n.z^/fp.u  -  p.)] 

-[njCPi  -  Pi)Zi(l  -  2P.)]/[P.(1  -  P.)]2 

The  essential  trick  in  the  approximate  Newton  iteration  (see  Fisher 
[5])  consists  of  ignoring  small  terms  in  the  second  derivative.  Thereby 
the  iterative  procedure  is,  admittedly,  slower  than  the  complete  Newton 
procedure.  However,  since  no  approximation  is  employed  in  the  first  deri¬ 
vatives,  the  final  result  will  still  be  exact  when  the  process  has  con¬ 
verged.  In  the  above  example,  (pi  -  P^)  is  a  relatively  small  term  and 
will  be  disregarded  in  the  second  derivative.  Also, 

3(31og  L/3u)/3Z. 

involves  such  a  difference  (p^  -  P^)  and  hence  will  be  disregarded. 

With  this  simplification  one  obtains  the  approximate  value 

k  k 

3 2 log  L/3y2  »  -  l  n.{Z2/[P. (1  -  P .)]}  *  -  [  n  W  .  (2.9a) 

i«l  i=l 

(2.9b) 

(2.9c) 


Similarly, 


32log  L/3B2 


r  —  2 

-A  n.W. (x.  -  x) 

i.i 111 


and 


32log  L/3u36  *  -  l  n.W. (x.  -  x) 
i-1  1  1  1 


■  0  if  x  is  chosen  as 
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k  k 

l  ;I.W.  X.  /  y  n.W.. 
.**,  ill  . L ,  11 


With  this  simplification,  the  Hessian  matrix  of  second  derivatives  is 

diagonal.  Hence  corrections  can  be  made  on  each  of  the  two  unknowns 

separately,  by  substitution  into  the  customary  Newton  formula. 

k  *  k 

Pi  =  pn  ♦  \  n.W.[(p.  -  P.)/Z.]/  n.W. 

1  0  >.  l  il'ri  i"  iJ'.L,  i  i 


k  k 

6n  +  y  n.W. [(p  -  P. )/Z.](x.  -  x)/  y  n.W.(x.  -  x)^ 
0  i  l i  l  iJV  i  i-\  111 


If  we  estimated  ard  8^  by  weighted  regression  on  Y^  given  x^ 
with  the  weights  nW,  we  could  combine  the  pQ  and  the  sum  into  a  single 
sum 

"i  ■  *  (pi  -  pi)/zi!'xi  si/JjW'i  -  *>2- 

and,  analogously 


y  n.W. [Y.  ♦  (p. 

l  il  l  v*i 

i=l 


P. )/Z. ] (x.  -  x)/  y  n.W. (x. 
1  iJ  l  *  ■  x  i  i 
i=l 


But  these  are  estimates  of  and  Bj  of  a  weighted  regression 
analysis,  between  the  x^  and  the  quantities  [Y^  +  (p^  - 
The  latter  are  called  "working"  quantities.  Thus,  the  approximate 
Newton  iteration  employs  the  same  technique  as  the  iterated  (weighted 
regression)  method,  with  y^  replaced  by  the  working  quantities, 
y.  +  (p.  -  P.)/Z.. 

7i  v*i  i 
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The  c.d.f.  most  frequently  employed  in  the  analysis  of  such  data  is 
the  Normal  distribution  function,  often  with  the  mean  placed  at  5  to 
avoid  negative  numbers.  This  is  called  the  Probit  transformation.  In 
detail,  the  following  steps  are  required: 

Data  for  the  analysis  are,  for  each  of  k  batches,  x^,  the  stimulus 
or  dosage  applied;  n^,  the  number  of  subjects  tested  at  each  dose;  and 
r^,  the  nuntoer  of  subjects  who  responded  at  each  dose. 

The  first  cycle  proceeds  as  follows: 

(a)  Compute  the  proportion  of  subjects  responding  to  various  levels  of 
a  stimulus  p.  =  r./n.  where  i  =  l,2,...k  levels  of  a  stimulus. 

(b)  Obtain  y^,  the  empirical  probit  for  each  dosage  level;  y^^  =<t_1(p.), 
where  #-1(p^)  is  the  inverse  normal  distribution  (e.g.,  <J>~ 1 C - 9^75) 

=  1.96  or  6.96  if  5  is  added). 

(c)  Initialize  the  working  probit  by  setting  y?  *  y^  and  the  weighting 
coefficient  for  the  probit  analysis  *  by  using  =  n^. 

(d)  From  the  weighted  regression  scheme,  obtain  a  and  b,  estimates  of 
a  and  0  in  the  equation  y^  *  a  +  0x^. 


b  =  s£)/s»  “here 


*  z  »,vi  ■  (( i  Vi,(  1  vi))/ 1  “i 

y  i-1  1  1  1  i-1  1  1  i=l  11  i=l  1 


,  k  .  k  ,  k 

s(w>  ■  y  w.x.  -  (  y  w.x.)  /  y  w.. 

**  ini  1  1  ini  1  ^  itl  i 


and 


a 


bx  where 


=  y 

k  k 

y  *  I  w.y./  [w.  and 
i=l  1  1  i=l 


x 


k 


w.x. 
1  l 


/  £w. . 
i-11 


(Note  that  w^  is  nJih  here). 

(e)  Compute  values  of  expected  probit  using  the  initial  estimates  of 
and  b:  Y.  =  a  +  b  x.  where  i  =  l,2,...k. 

(f)  Calculate  the  predicted  proportion  of  success  at  each  dose. 
Here,  =  *(Yi)[or,  possibly,  Y^  -  5],  is  the  cumulative  normal 
distribution  value  corresponding  to  Y^. 

(g)  Obtain  Z^  *  fCY.),  the  ordinate  under  the  normal  distribution 
curve,  for  each  Y^. 

(h)  Corresponding  to  each  value  of  Y,  compute  the  following: 

W.  =  Z?/[P.(1  -  P.)]  and  y*.  *  Y.  +  (p.  -  P.)/Z.. 
l  r  1  iJi  'i  i  ri  i  i 

(i)  Return  to  step  (d)  and  begin  a  new  iteration  with  a  ■  a  and  b  =  8 

until  the  relative  difference  between  two  consecutive  values  of  a 

-4 

or  8  are  small.  Relative  error  of  10  ,  or  5  iterations,  are 

employed  in  this  program. 

Other  models  may  be  substituted  in  this  analysis.  The  following 
equations  for  P,  y,  and  Z  should  be  substituted  in  steps  (b),  (f),  and 
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The  logistic  function,  dn/dt  =  Xn(c  -  n),  where  c  is  the  ceiling  and 

n  is  some  population  size,  is  widely  used  to  represent  such  phenomena 

Y  Y  -Y 

as  population  growth.  P  *»  e  /(I  ♦  e  )  =  1/(1  +  e  )  gives  the  metametric 

equation  and  Y'=  log[P/(l  -  P)].  So  that  Z  =  dP/dY  =  eY/(l  +  eY)2. 

The  log-log  or  "Gumpertz"  function  is  used  when  reproduction  inhibits 

-e-Y 

itself.  The  metametric  equation  is  P  =  e  and  Y  =  -log(-logP). 

Also  Z  -  dP/dY  >  -P  log  P. 

For  the  Arc  Sine  function,  P  =  fs(l  +  sin  Y) ,  Y  «  Arc  sin  (2P  -  1), 
and  Z  =  Jscos  Y. 

-Yr  1/r 

Finally,  for  the  Weibull  function,  P  =  1  -  e  ,  Y  =  [-log  (1  -  P)]  , 


CHAPTER  III 

THE  INTERACTIVE  PROGRAM 


Programmed 
Function 
Keyboard 

The  diagram  above  displays  the  position  of  the  IBM  2250  screen,  the 
programmed  function  keyboard,  and  the  alphameric  keyboard. 

The  first  message  on  the  screen,  given  by  the  GMS  monitor,  states 
that  the  iser  must  depress  any  one  of  the  programmed  function  keys. 

The  second  message  lists  several  options.  The  user  should  first 
type  $NAMES,  then  depress  the  "Altn  Coding"  key  and  while  it  is  held 
down,  depress  the-  "5"  key  on  the  alphameric  keyboard.  The  list  of 
•programs  should  be  examined  carefully  to  ascertain  that  the  QUANTAL 
program  is  listed.  Next,  one  should  type  $LINK  QUANTAL,  followed  by 
the  "Altn  Coding"  and  "5"  key. 

Once  the  QUANTAL  program  has  been  successfully  obtained,  a  message 
will  be  displayed  and  the  user,  if  entering  for  the  first  time,  should 
depress  programmed  function  key  1.  If  the  user  is  entering  the  program 
following  a  departure  to  the  CALCG  program  (see  below) ,  he  should  depress 
programmed  function  key  2  to  return  to  his  point  of  departure. 
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Next,  the  general  instructions  of  the  program  will  be  listed. 

These  instructions  are  important.  They  state: 

1.  At  any  point  in  the  program  you  may  restart  by  depressing 
programmed  function  key  30  or  terminate  by  depressing  programmed 
function  key  31. 

2.  Your  answers  should  be  typed  from  the  typewriter  keyboard  direct¬ 
ly  in  front  of  you  (unless  you  are  asked  differently). 

3.  All  numbers  are  to  be  typed  as  real  numbers  --  decimal  point 
must  follow  the  number. 

4.  After  each  statement  you  should  cause  an  "End  of  Block  (EOB)" 
signal,  by  first  depressing  the  "Altn  Coding"  key  and,  while  it 
is  held  down,  depressing  the  "5"  key. 

5.  If  there  is  more  than  one  blank  to  be  filled  in  a  line,  depress 
the  "JUMP"  key  located  on  the  left  side  of  the  typewriter  key¬ 
board  to  proceed  from  blank  to  blank. 

The  programmed  function  key  1,  employed  throughout  the  program  to 
relate  the  normal  flow  of  the  program,  should  be  depressed  when  the  user 
is  ready  to  continue. 

Now  the  user  will  be  asked  to  enter  his  first  piece  of  data:  the 
number  of  classes  or  batches  on  v/hich  he  wishes  to  perform  the  analysis. 
This  number  must  be  greater  than  zero  and  less  than  one  hundred.  The 
number  of  classes  or  batches  will  ordinarily  range  between  5  and  25. 

With  fewer  than  4  classes  the  results  are  probably  meaningless;  more  than 
25  classes  would  be  meaningful  only  if  the  total  sample  were  very  large. 

In  the  event  that  a  number  is  entered  which  is  not  in  the  required  range, 
an  error  message  will  be  displayed  and  the  user  must  re-enter  a  correct 


13 


number.  The  user  must  note  that  each  number  must  be  entered  as  a  real 
number  (with  a  decimal  point  actuall/  typed,  even  if  the  number  is  an 
integer) . 

The  next  sequence  of  statements  will  request  information  on  the 
concomitant  variable  or  x-dosage  which  was  administered  to  each  batch, 
the  number  of  individuals  or  subjects  in  each  batch,  and  the  number  of 
subjects  in  each  batch  which  responded.  These  three  values  will  be  filled 
into  the  blanks  of  the  sentence:  In  batch  xx  where  the  concomitant 

variable  (dosage)  _  was  given  to  _  subjects,  there  were 

_ subjects  that  responded  (successes).  The  user  may  proceed  from 

blank  to  blank  by  depressing  the  "JUMP"  key.  The  "End  of  Block"  should 
be  caused  only  when  the  entire  sentence  has  been  completed.  In  the 
first  blank,  the  x-dosage  (the  mean  x-score)  for  each  batch  may  be 
entered  as  a  real  number  in  F  or  E  (scientific  notation)  field  specifica¬ 
tion.  The  width  of  the  field  is  14,  thus  data  may  be  entered  with  maxi- 

mums  F14.7  or  E14.7.  If  E  format  is  used,  the  exponent  (decimal  point 

-8 

shift)  must  have  two  digits.  For  example,  1.62  *  10  would  be  entered 
as  1.62E-08.  The  x-dosage  may  be  either  negative  or  positive.  However, 
since  one  of  the  important  aspects  of  this  program  is  to  permit  trans¬ 
formations  on  x  -  in  order  to  assure  better  fit  -  the  dosages  should, 
as  a  rule,  be  positive  (log  ahd  square-root  are  frequently  used  trans¬ 
formations)  . 

In  the  second  blank,  the  number  of  subjects  or  individuals  in  each 
batch  is  requested.  It  must  be  greater  than  two  and  must  be  entered  as 


a  real  number. 
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The  number  of  subjects  which  responded  in  each  class  occupies  the 
third  blank.  It  must  be  zero  or  greater  and  must  not  exceed  the  number 
of  individuals  in  each  batch. 

If  any  of  these  requirements  are  not  fulfilled,  an  error  message 
will  be  displayed  and  the  user  must  complete  the  statement  again. 

When  the  last  statement  of  input  has  been  successfully  entered, 
the  data  will  be  displayed  in  an  array.  The  user  should  examine  his 
data  carefully  to  ascertain  that  all  data  have  been  entered  correctly. 

The  user  will  then  be  given  the  opportunity  to  change  data  of  any  of  the 
batches,  to  eliminate  a  batch,  to  add  a  new  batch,  or  to  continue  with 
the  analysis. 

If  the  user  wishes  to  correct  any  of  the  batches  (more  than  one 
batch  may  be  altered),  he  should  depress  programmed  function  key  2.  At 
this  time  a  message  will  appear  that  asks  which  batch  the  user  desires 
to  alter.  If  the  number  typed  in  by  the  user  is  not  in  the  correct  range, 
an  error  message  will  appear  and  the  user  must  state  the  correct  batch 
number.  When  this  number  has  been  entered  satisfactorily,  a  statement 
to  obtain  inpui.  will  be  displayed  on  the  screen.  The  user  must  complete 
this  statement  just  as  when  he  entered  it  the  first  time.  However,  he 
may  eliminate  this  batch  by  simply  refusing  to  fill  in  any  of  the  blanks 
(responding  EOB  at  once). 

If  the  user  wishes  to  add  another  batch  to  his  data,  he  should 
depress  programmed  function  key  5.  It  will  then  be  necessary  to  complete 
the  input  statement  according  to  the  requirements  pertaining  to  it. 

After  an  addition  or  alteration  has  been  made  by  the  user,  the 
entire  data  will  again  be  displayed  --  including  corrections  and  addi- 
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tions  --  and  the  user  has  the  option  once  again  to  add  or  alter  a  batch 
or  to  continue. 

When  the  data  completely  satisfy  the  user,  he  should  depress  pro¬ 
grammed  function  key  1.  At  this  time  the  analysis  will  be  initiated 
and  the  user  must  decide  which  transformation  he  desires  for  his  data. 

He  may  choose  the  PROBIT,  LOGIT,  LOG-LOG,  ARC  SIN,  or  WEIBULL  by  depres¬ 
sing  programmed  function  key  6,  7,  8,  9,  or  10  respectively.  In  the 
WEIBULL  case  he  will  also  need  to  supply  the  exponent  by  depressing  the 
appropriate  key. 

At  this  stage  in  QUANTAL,  an  alarm  will  be  sounded  and  the  screen 
will  be  blank  for  a  few  seconds.  It  is  somewhat  unpredictable  how  long 
this  will  take.  The  GMS  monitor  is  time-shared.  If  the  processor  queue 
is  very  long,  a  minute  or  two  may  pass.  If  there  is  little  batch  activi¬ 
ty,  a  few  seconds  may  be  all  that  is  needed. 

The  first  display  of  the  output  section  will  be  the  graph  of  the 
observed  and  predicted  proportion  of  successes  for  each  class.  The 
x-dosages  have  been  arranged  in  ascending  order  and  are  depicted  on 
the  horizontal  axis.  The  batch  numbers  are  written  directly  below  this 
axis  and  the  minimum  and  maximum  x-dosage  levels  are  written  below  the 
respective  batch  numbers.  The  vertical  axis  has  the  proportion  of  suc¬ 
cesses.  The  lined  histogram. is  the  predicted  proportion  of  successes 
and  the  observed  proportions  of  success  are  denoted  by  asterisks.  The 
user  should  depress  programmed  function  key  1  when  he  wishes  to  continue. 

The  next  display  (or  series  of  displays)  contains  numerical  answers 
(ten  batches  per  display).  On  the  final  "page,"  the  user  may  return  to 
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"page  1"  by  depressing  programmed  function  key  4.  Displayed  for  each  of 
the  batches  is  the  size  of  the  batch,  P(OBS)  —  the  observed  proportion 
of  successes,  P(PRED)  —  the  predicted  proportion  of  successes,  Y(PRED) 

—  the  predicted  y-scores,  Y(OBS)  —  the  observed  y-scores,  Y(WORK)  — 
the  "working"  quantities,  the  weights,  X-DOSAGE  --  the  observed  x-scores, 
and  Y-ERROR  --  the  error  of  prediction  in  the  y's. 

The  next  display  in  the  output  section  is  an  analysis  of  variance 
table.  The  sum-of-squares,  degrees  of  freedom,  and  mean-squares  are 
given  for  the  regression  and  error.  Also  the  calculated  F  =  mean-squares 
regression/mean-squares  error  is  displayed.  The  user  then  has  the  option 
to  leave  the  QUANTAL  in  order  to  call  the  CALCG  program,  e.g.,  to  find 
whether  or  not  his  F  value  is  significant. 

Upon  returning  to  the  QUANTAL  program,  if  the  user  chose  to  call 
the  CALCG  program,  he  is  reminded  to  return  to  this  point  by  depressing 
programmed  function  key  2. 

If  the  user  does  not  wish  to  leave  the  QUANTAL  program,  he  should 
depress  programmed  function  key  1  when  he  wishes  to  continue. 

A 

Next  there  is  a  graph  of  the  regression  line  Y  =  &  +  Bx  where 
d  and  B  have  been  estimated  by  the  appropriate  quantal  analysis,  and 
x  is  the  dosage.  The  horizontal  axis  is  the  dosage-axis  and  the  verti¬ 
cal  axis  is  the  y-axis.  The  observed  y-values  are  denoted  on  the  graph 
by  asterisks.  Programmed  function  key  1  is  depressed  to  continue. 

In  the  next  section  the  user  may  alter  his  x-dosage  by  one  of  the 
following  functions:  sine,  cosine,  tangent,  arc  sine,  arc  cosine, 
arc  tangent,  exponential,  natural  logarithm,  common  logarithm,  square 
root,  absolute  value  and  x-squared.  To  use  one  of  the  fun„t?ons,  he 


depresses  programmed  function  key  3  and  then  types  in  the  number  corres¬ 
ponding  to  the  function  desired.  If  a  function  is  chosen  which  is  in¬ 
compatible  with  the  data,  e.g.,  square-root  and  negative  x-dosages,  an 
error  message 'is  displayed  and  the  user  must  give  correct  answers. 

If  a  function  is  chosen  to  alter  the  x-dosages,  then  the  next  sec¬ 
tion  will  be  a  display  of  the  final  (regression  line)  section  with  the 
x-scale  altered.  In  fact  the  line  itself  will  not  be  displayed  --  only 
the  asterisks.  If  the  user  likes  this  new  plot  he  may  perform  reanalysis 
by  continuing  (depressing  programmed  function  key  1).  If  he  does  not 
like  it  he  will  be  given  a.  chance,  in  the  next  display,  to  call  the 
inverse  of  the  first  function  and  then,  after  yet  another  return,  try  a 
different  function. 

At  the  stage  where  the  user  is  asked  to  supply  the  transformation 
he  has  an  option  to  depress  programmed  function  key  1  and  thus  obtain 
a  new  analysis  (starting  with  the  Summary)  of  the  data  in  accordance 
with  the  last  chosen  transformation. 


18 


ILLUSTRATION 

The  QUANTAL  program  is  called  by  the  user  and  the  instructions  as 
outlined  in  the  User's  Guide  are  followed. 

The  data  for  the  analysis  are  typed  in  by  the  user  by  his  completing 
the  sentences  as  shown  in  figure  3.1. 

When  all  data  have  been  entered,  a  tabular  display  of  the  data  is 
exhibited  (figure  3.2). 

In  this  illustration,  the  probit  transformation  has  been  chosen  by 
the  user's  depressing  programmed  function  key  6. 

After  a  short  delay,  for  the  calculations,  there  appears,  on  the 
screen,  a  step  function  indicating  the  estimated  cumulative  distribution 
function  (proportion  vs.  dosage).  The  asterisks  show  the  observed  values. 
If  the  fit  were  perfect,  they  would  appear  coincident  with  the  upper  cor¬ 
ner  of  each  step.  Figure  3.3  shows  this  display. 

Next  there  appears  a  detailed  numerical  display  (figure  3.4). 

Notice  here  that  the  greatest  y-error  is  -.25310  in  batch  4. 

The  next  display  (figure  3.5)  is  an  analysis  of  variance  table. 

Next,  the  graph  of  the  straight  line  Y  ■  a  +  6*  (as  shown  in  figure 
3.6)  is  displayed.  It  will  be  noted  here  that  the  analysis  produces 
a  very  good  fit  for  all  but  one  of  the  data  points. 

In  the  second  stage  of  this  presentation,  the  user  does  not  alter 
the  x-dosages  but  chooses  the  logit  transformation. 

Figures  3.7  and  3.8  display  the  numerical  results  and  the  analysis 

of  variance  for  the  logistic  transformation.  As  a  comparison  of  the 

2 

closeness  of  fit  one  could  calculate  R  ■  sum-of-squares  regression/ 
sum-of-squares  Total  which  is  .9881  for  the  Probit  and  .9866  for  the 
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logistic.  This  is  an  example  of  the  (usual)  closeness  of  the  normal 
and  logistic  transformations.  Note  that  no  other  base  for  comparison 
would  be  possible,  since  the  dispersion  of  probity  and  logits  are  differ 
ent  (.975  corresponds  to  a  probit  of  1.96  (6.96)  but  to  a  logit  of  3.66) 


FIGURE  3.1 

DATA  INPUT  STATEMENTS 
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FIGURE  3.2 

TABULAR  DISPLAY  OF  DATA 
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THE  ESTIMATED  CUMULATIVE  DISTRIBUTION 
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NUMERICAL  RESULTS  OF  THE  ANALYSIS 
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ANALYSIS  OF  VARIANCE  TABLE 
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NUMERICAL  RESULTS  OF  THE  ANALYSIS 
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ANALYSIS  OF  VARIANCE  TABLE 


LOGIT  ANALYSIS 
ANALYSIS  OF  VARIANCE  TABLE 


NOT  REPRODUCIBLE 


CHAPTER  IV 

PROGRAM  DOCUMENTATION 

The  QUANTAL  program  has  been  written  in  FORTRAN  IV.  The  listing  is 
presented  in  Appendix  A.  Inasmuch  as  execution  of  the  program  must  use 
no  more  than  64  K  bytes  of  core  (16  K  words),  an  overlay  structure  is 
required.  The  exposition  (figure  4.1)  shows  the  job  control  cards  for 
this  overlay. 

As  may  be  seen  in  figure  4.2,  the  overlay  structure  for  the  QUANTAL 
program  is  relatively  simple.  The  controlling  segment,  or  root  segment 
contains  the  MAIN  program  and  several  basic  GRAF  (Graphics  Addition  to 
FORTRAN)  routines.  This  segment  remains  in  storage  throughout  execution 
of  the  QUANTAL  program  and  branches  off  into  two  separate  segments.  One 
of  the  segments  contains  the  FORTRAN  routines  that  are  necessary  for  the 
mathematical  calculations  of  the  program  and  also  the  subroutines  CALCG 
and  FNS.  The  other  branch  of  QUANTAL  contains  the  GRAF  routines  and  the 
subroutines  which,  through  use  of  these  routines,  obtain  the  data  from 
the  user  and  display  the  results  of  the  analysis.  The  three  subroutines, 
OUTPUT,  CHNGTR,  and  GRACDF,  require  the  GRAF  routines,  CORDCALL,  LINE$$, 
SIZE,  $CORD$,  and  POINTS,  to  present  various  graphical  displays.  Only 
one  subroutine,  with  the  necessary  FORTRAN  or  GRAF  routines  provided  by 
the  structure  resides  in  storage  with  the  root  segment  at  any  one  time 
and  thus  the  total  size  of  the  program  i»  appropriately  restricted. 
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THE  MAIN  PROGRAM 

The  reader  is  referred  to  reference  [7]  for  an  explanation  of  the 
meaning  of  each  of  the  GRAF  (Graphics  Addition  to  FORTRAN)  subroutines. 

The  main  program  of  the  QUANTAL  program  is  designed  to  call  various 
subroutines,  in  an  order  primarily  designated  by  the  user.  The  flow  of 
the  program  is  determined  by  the  value  of  KEY,  a  fixed  point  value  which 
is  stored  in  COMMON. 


VALUE  OF  KEY  COMMENT 


1 

2 

3 

4 

5 

6 
7 

10 

50 


Restart  the  program 
Proceed  to  display  data 
*  Proceed  to  choose  transformation 
Proceed  to  change  function 
Proceed  to  display  data 
Proceed  to  alter  x-dosage 
Same  as  3 

Alter  one  class  of  data 
Terminate  the  program 


The  COMMON  set 


NDET: 

An  array  required  by  GRAF 

X'(100) : 

Reserved  for  x-dosage 

N(100) : 

Number  of  subjects  per  batch 

NR(100) : 

Number  of  successes  per  batch 

DV1  to  DVB 

inclusive:  DV  (display-variable)  names 

NCL: 

Number  of  batches 

KEY: 

Control  of  program  flow  (see  above) 

KK: 

Designation  for  batch  to  be  changed 

IFUN: 

Indicates  function  on  x-dosage 

ITRAN: 

Indicates  tranformation  on  P  (Probit 

IREC: 

Index  for  random  access  I/O. 

Subroutine  inputa 
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Purpose:  The  purpose  of  subroutine  INPUTA  is  to  obtain  the  number  of 
classes  in  the  analysis. 


Procedure: 

FORTRAN 

STATEMENT  STATEMENT 

NUMBER  NUMBER 

3 


4 

6-8 

9-12 


13  402 

14 


COMMENT 

TTie  Graphics  Data  Control  Block  (GDCB)  is  opened 
and  the  display  variables  DV1,  DV2,  DV3,  DV4, 
and  DVE  are  established,  thereby  initializing 
the  2250  for  I/O. 

The  lights  of  programmed  function  keys  1-11  and 
30  and  31  are  turned  on. 

The  scope  and  display  variables  are  cleared  from 
previous  commands. 

The  message  in  FORMAT  400  is  plotted  on  the  screen 
to  enable  the  user,  if  returning  from  the  CALCG 
program  to  proceed  directly  to  his  point  of 
departure  in  the  QUANTAL  program. 

Procedure  is  detained  until  the  user  responds. 

If  programmed  function  key  2  has  been  depressed, 
(NDET(4)  equals  2),  indicating  a  return  from 
the  CALCG  program,  proceed  to  FORTRAN  statement 


15 

16 


401. 

If  programmed  function  key  1  has  been  depressed, 
proceed  to  FORTRAN  statement  66. 

If  neither  of  these  programmed  function  keys 
has  been  depressed,  return  to  statement  402 
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SUBROUTINE  INPUTA  (Continued) 


FORTRAN 

STATEMENT  STATEMENT  COMMENT 

NUMBER  NUMBER 


17  66 

18 

19-25 

26  201 

27 

28 

29 

30 

31-34  25 


(await  proper  reply). 

The  scope  is  cleared  of  all  existing  displays. 
Display  variable  DV1  is  reinitialized. 

The  messages  in  FORMAT  statements  100  and  1000, 
the  general  instructions  of  the  program,  are 
displayed  on  the  screen. 

Procedure  is  detained  until  user  has  finished 
reading  the  display. 

If  programmed  function  key  31  has  been  depressed 
(NDET(4)  equals  31),  indicating  the  user  wishes 
to  terminate  the  program,  proceed  to  FORTRAN 
statement  750  (go  out) . 

If  programmed  function  key  30  has  been  depressed 
(NDFT(4)  equals  30),  indicating  the  user  wishes 
to  restart  the  program,  proceed  to  FORTRAN 
statement  67. 

If  an  interrupt  has  been  given  by  the  user  from 
any -other  device  except  the  programmed  function 
keyboard  (I  is  not  equal  to  1),  return  to  detain 
statement  201. 

The  scope  is  cleared  of  all  existing  displays. 
The  message  in  FORMAT  statement  101,  asking 
the  user  to  type  in  the  number  of  batches  in 
the  analysis,  is  plotted  on  the  screen. 
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FORTRAN 

SUBROUTINE  INPUTA  (Continued) 

STATEMENT 

STATEMENT 

COMMENT 

NUMBER 

NUMBER 

35 

• 

A  cursor  is  placed  in  the  first  unprotected 

• 

position  of  FORMAT  statement  101  so  that  the 

user  can  type  in  the  number  of  batches. 

36 

200 

Procedure  is  detained  until  user  types  in 

number  of  batches. 

37 

If  programmed  function  key  31  has  been  depressed, 

indicating  the  user  wishes  to  terminate  the  pro¬ 
gram,  proceed  to  FORTRAN  statement  750. 

38 

If  programmed  function  key  30  has  been  depressed. 

indicating  the  user  wishes  to  restart  the  pro¬ 

gram,  proceed  to  FORTRAN  statement  67. 

39 

69 

If  an  interrupt  has  been  given  by  the  user  from 

- 

any  device  except  the  alphameric  keyboard  (I  is 

not  equal  to  2)  return  to  detain  statement  200. 

40-43 

XNUM  is  NCL,  floated. 

44 

The  pointers  for  dummy  unit  4  are  reset. 

45-47 

The  cursor  is  removed  from  the  screen. 

48 

The  integer  variable  NCL  is  equal  to  XNUM. 

49 

If  the  number  of  classes  (NCL)  is  less  than  zero 

or  greater  than  100,  proceed  to  FORTRAN  state¬ 
ment  30  where  an  error  message  is  displayed. 

50 

50 

When  the  statement  has  been  completed  correctly 

all  error  messages  are  erased. 

SUBROUTINE  INPUTA  (Continued) 


41 


FORTRAN 

STATEMENT  STATEMENT  COMMENT 

NUMBER  NUMBER 

52-55  30  The  error  message  in  FORMAT  statement  114  is 

displayed  on  the  screen. 

56-57  The  incorrect  statement  101  is  removed  from  the 

Buffer  Table  and  display  variable  DV2  is  re¬ 
initialized. 

58  Display  variable  DVE  is  re-initialized. 

60  401  KEY  is  set  equal  to  25. 

61  The  screen  is  cleared  of  all  existing  displays. 

64  21  Display  variables  DV1,  DV2,  DV3,  DV4,  and  DVE 

are  re-initialized. 


SUBROUTINE  DATA 


Purpose: 

Procedure 

STATEMENT 

NUMBER 

3 

4 

5 

6 

7 

8 

9-12 

13-19 

20-23 


The  purpose  of  this  subroutine  is  to  obtain  tha  data  from  the 
user  by  his  completing  the  sentence:  In  batch  XX  where  the 

concomitant  variable  (dosage)  _  was  given  to  _ 

subjects,  there  were  _  subjects  that  responded  (successes). 


FORTRAN 

STATEMENT  COMMENT 
NUMBER 

The  display  variables  DV3,  DV4,  and  DVE  are 
established. 

The  scope  is  cleared  of  all  existing  displays. 
If  KEY  is  equal  to  10,  indicating  data  is  to 
be  obtained  from  only  one  batch  (when  changes 
are  made),  proceed  to  FORTRAN  statement  75. 

NL,  the  number  of  times  FORMAT  statement  103 
is  to  be  written  is  initialized  equal  to  NCL. 
JEY  is  a  flag  =  0  for  all  but  the  last  page, 

*  1  for  the  last  page. 

NI,  the  increment  of  loop  500  is  initialized. 

4  The  message  in  FORMAT  102  to  obtain  data  from 

a  batch  is  plotted  on  the  screen. 

Only  eight  statements  are  to  appear  on  the 
screen  at  a  time.  NN  is  the  delimiter  of  loop 
500. 

The  statement  to  obtain  data  from  the  user  is 
plotted  on  the  screen. 
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SUBROUTINE  DATA  (Continued) 


FORTRAN 

STATEMENT  STATEMENT 

NUMBER  NUMBER 

24 


25  201 


26 


27 


•  28 


29-32 


33 

34 


COMMENT 

A  cursor  is  placed  in  the  first  unprotected 
position  of  FORMAT  statement  103  so  that  the 
user  can  type  in  data. 

Procedure  is  detained  until  the  user  has  entered 
data  for  a  batch. 

If  programmed  function  key  (PFK)  31-  has  been 

depressed  (NDET(4)  equals  31),  indicating  the 

user  wishes  to  terminate  the  program,  proceed 

to  FORTRAN  statement  750  (exit).  j 

If  PFK  30  has  been  depressed  (NDET(4)  equals  > 

1 

30),  indicating  the  user  wishes  to  begin  the  ! 

program  again,  proceed  to  FORTRAN  statement  67. 

If  an  interrupt  has  been  given  by  the  user  from  j 

any  device  except  the  alphameric  keyboard  (I  is  j 

1 

not  equal  to  2)  return  to  detain  statement  201. 

X(KK),  the  x-dosage;  and  XX  and  XY,  the  real 
values  for  the  number  of  individuals  in  batch 
KK-  and  the  number  of  successes  in  batcn  KK, 
are  read  according  to  FORMAT  301. 

The  pointers  for  dummy  unit  4  are  reset. 

N(KK)  is  defined  equal  to  XX.  ! 

I 

NR(KK)  is  defined  equal  to  XY. 


35 
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Si  M  Oim  :  DAT/  .  tinued) 

FORTRAN 

STATEMENT  STATEMENT  COMMF.tfl 

NUMBER  NUMBER 

37  -22  If  the  number  of  individuals  in  each  batch  is 

less  than  3,  proceed  to  FORTRAN  statement  11  to 
write  an  error  message. 

38  12  If  the  number  of  individuals  who  responded  is 

negative,  proceed  to  FORTRAN  statement  13  to 
write  an  error  message.  If  the  number  of 
individuals  is  positive,  proceed  to  next  state¬ 
ment. 

39  15  If  the  number  of  individuals  who  responded  is 

greater  than  the  number  of  individuals  in  that 
batch,  proceed  to  FORTRAN  statement  45  to  write 
an  error  message. 

■41-43  20-21  The  beam  is  placed  in  the  upper  left-hand  comer 

of  the  screen. 

45-48  11  The  error  message  in  FORMAT  statement  117  is 

displayed  on  the  screen. 

49-50  The  incorrect  statement  103  is  removed  from 

the  Buffer  Table  and  display  variables  DV4  and 
DVE  are  re-initialized. 

52-55  13  The  error  message  in  FORMAT  statement  118  is 

plotted  on  the  screen, 

56-57  The  incorrect  statement  103  is  removed  from 

the  Buffer  Table  and  display  variables  DV4  and 


DVE  are  re-initialized. 
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SUBROUTINE  DATA  (Continued) 


STATEMENT 

FORTRAN 

STATEMENT 

COMMENT 

NUMBER 

NUMBER 

59-62 

45 

The  error  message  in  FORMAT  statement  120  is 

• 

displayed  on  the  screen. 

63 

• 

Incorrect  statement  103  is  removed  from  the  BT. 

64 

The  cursor  is  removed  from  the  screen. 

65 

Display  variables  DV4  and  DVE  are  re-initialized. 

67-69 

501 

The  cursor  is  removed  from  the  screen. 

70 

Display  variable  DV4  is  re-initialized. 

71 

500 

End  of  loop  500  -  All  error  messages  are  erased 

72 

from  the  screen. 

The  scope  is  cleared  of  all  existing  displays. 

73-74 

NI  and  NL  -  the  increment  and  decrement  of  loop 

"75 

500  are  up-dated. 

If  JEY  is  equal  to  1,  indicating  all  data  has 

76-77 

75 

been  obtained,  exit;  otherwise,  proceed  to 

FORTRAN  statement  4. 

When  subroutine  DATA  has  been  called  with  KEY=10, 

78 

the  increment  and  delimiter  of  loop  500  are  set 

to  the  desired  batch  KK. 

JEY  is  set  equal  to  1. 

79-82 

FORMAT  statement  192  is  displayed  on  the  screen. 

84 

750 

KEY  is  set  equal  to  50. 

86 

67 

KEY  is  set  equal  to  1. 

87 

99 

Display  variables  DV3,  DV4,  and  DVE  are  rc-ini- 

tialized. 
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Purpose; 

Procedure: 

STATEMENT 

NUMBER 

3 

4 

5-7 

8 

9 


10-13 

14 

15 


SUBROUTINE  INPUTB 

The  purpose  of  this  subroutine  is  to  display  the  summary  of 
the  data  obtained  in  subroutine  DATA  and  to  allow  the  user  to 
alter  any  data  or  to  add  or  remove  one  or  more  classes  to  the 
analysis. 

FORTRAN 

STATEMENT  COMMENT 
NUMBER 

The  display  variables  DV5,  DV6,  DVT,  DVD,  DVA, 

•  and  DVE  are  established. 

502  The  scope  is  cleared  of  all  existing  displays. 

The  title  message,  FORMAT  statement  104,  is 
written  onto  display  variable  DV5. 

K,  the  number  of  batches  to  be  displayed  on  the 
screen  at  one  time,  is  initialized  equal  to  0. 

Loop  21  -  One  class  of  individuals  at  a  time, 

and  the  corresponding  x-value,  N-value,  and 
NR-value  are  to  be  plotted  on  the  screen  at  a 
time. 

The  number  of  the  class,  the  x-value,  N-value, 
and  NR-value  are  plotted  on  the  screen  according 
to  FORMAT  statement  105. 

Display  variables  DV5,  and  DV7  are  re-initialized. 
K  is  incremented  by  one. 

A  statement  to  insure  that  no  more  than  25  rows 
appear  on  a  single  display. 


15 
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STATEMENT 

FORTRAN 

STATEMENT 

SUBROUTINE  INPUTB  (Continued) 

COMMENT 

NUMBER 

NUMBER 

17 

. 

The  screen  is  cleared  of  all  existing  displays. 

18 

K  is  re-initialized  equal  to  0. 

19-20 

. 

The  title  message,  FORMAT  statement  104,  is 

22-25 

written  onto  display  variable  DV5. 

A  message  which  gives  the  user  the  ability  to 

correct  any  data  or  to  add  another  class  to  his 

26 

202 

analysis  is  displayed  on  the  screen. 

Procedure  is  detained  until  user  decides  whether 

27 

or  not  he  wishes  to  change  data. 

If  an  interrupt  has  been  given  by  the  user  from 

any  device  except  the  programmed  function  key¬ 
board  (J  is  not  equal  to  1)  return  to  detain 

statement  202. 

28 

7 

All  error  messages  are  cleared  from  the  screen. 

29 

If  programmed  function  key  (PFK)  31  has  been 

30 

6 

depressed  (NDET(4)  equals  31),  indicating  the 

user  wishes  to  terminate  the  program,  proceed 

to  FORTRAN  statement  710  (exit). 

If  PFK  30  has  been  depressed  (NDET(4)  equals  30), 

31 

indicating  tne  user  wishes  to  restart  the  program 

proceed  to  FORTRAN  statement  67. 

If  PFK  1  has  b«en  depressed  (NDET(4)  equals  1), 

indicating  the  user  does  not  wish  to  alter  any 
data,  proceed  to  FORTRAN  statement  76. 


SUBROUTINE  INPUTB  (Continued) 
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STATEMENT 

NUMBER 

32 

33 

34 

35 

36 

37 

38 

40 

41 

42-45 

46 


FORTRAN 

STATEMENT 

NUMBER 


300 


68 


COMMENT 

If  PFK  2  has  been  depressed  (NDET(4)  equals  2), 
indicating  the  user  wishes  to  make  a  correction 
to  his  data,  proceed  to  FORTRAN  statement  68. 

If  PFK  5  has  been  depressed  (NDET(4)  equals  5), 
indicating  the  user  wishes  to  add  another  batch 
to  his  analysis,  proceed  to  FORTRAN  statement  300. 
If  none  of  PFK  31,30,  1,  2,  or  5  have  been 
depressed,  return  to  detain  statement  202. 

The  number  of  classes  in  the  analysis,  NCL,  is 
increased  by  one  (  a  batch  is  to  be  added). 

The  scope  is  cleared  of  all  existing  displays. 

KK,  the  number  of  the  batch  to  be  cdded  is  set 
equal  to  NCL. 

The  variable  KEY  is  set  equal  to  10,  to  indicate 
subroutine  DATA  should  be  called  again. 

The  scope  is  cleared  of  all  existing  displays. 
Display  variable  DVA  is  re-initialized. 

The  message  in  FORMAT  statement  203,  which 
determines  in  which  batch  the  user  wishes  to 
change  the  data,  is  displayed  on  the  screen. 

A  cursor  is  placed  in  the  first  unprotected 
position  of  FORMAT  statement  203  so  that  the  user 
can  type  in  the  batch  number  to  be  corrected. 


63-66 


209 


50 


• 

SUBROUTINE  INPUTB  (Continued) 

FORTRAN 

• 

STATEMENT 

■  STATEMENT 

COMMENT 

'NUMBER 

NUMBER 

67 

.Incorrect  statement  203  is  cleared  from  the  screen. 

68 

Display  variable  DVE  is  re-initialized. 

70 

750 

KEY  is  equated  to  50. 

72 

76 

KEY  is  equated  to  7. 

74 

67 

KEY  is  equated  to  1. 

75 

74 

The  screen  is  cleared  of  all  existing  displays. 

76 

Display  variables  DV1,  DV2 ,  DV3,  DV4,  DV5,  i)V6, 

DV7,  DVD.  and  DVE  are  re- initialized. 

/ 

f 
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Purpose: 

Procedure: 

STATEMENT 

NUMBER 

3 

4-7 

8 

9 

10 

11 

12 


SUBROUTINE  CHNGTR 

Subroutine  CHNGTR  allows  the  user  a  choice  between  five  trans¬ 
formations:  the  Probit,  the  Logistic,  Log-Log,  Arc  Sin  and 
Weibul  1 . 

FORTRAN 

STATEMENT  COMMENT 

NUMBER 

Display  variable  DV1  is  established. 

FORMAT  statement  1  is  displayed  on  the  screen. 

2  Procedure  is  detained  until  the  user  has  chosen 

the  transformation. 

If  an  interrupt  has  been  given  by  the  user  from 
any  device  except  the  programmed  function  key¬ 
board  (J  is  not  equal  to  1)  return  to  detain 
statement  2. 

If  programmed  function  key  (PFK)  31  has  been 
depressed  (NDET(4)  equals  51),  indicating  the 
user  wishes  to  terminate  the  program,  proceed 
to  FORTRAN  statement  750  (exit). 

If  PFK  30  has  been  depressed  (NDET(4)  equals 
30),  indicating  the  user  wishes  to  restart  the 
program,  proceed  to  FORTRAh  statement  67. 

If  PFK  6  has  been  depressed  (NDET(s)  equals  6), 
indicating  the  user  desires  the  Probit  transfer 
mation,  ITRAN  is  set  equal  to  1. 


SUBROUTINE  CHNGTR  (Continued) 
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FORTRAN 

STATEMENT  STATEMENT 

NUMBER  NUMBER 


13 

• 

14 

15 

16 

17 

19 

750 

21 

67 

22 

99 

COMMENT 

If  PFK  7  has  been  depressed  (NDET(4)  equals  7), 
indicating  the  user  desires  the  Logit  transfor¬ 
mation,  ITRAN  is  set  equal  to  2. 

If  PFK  8  has  been  depressed  (NDET(4)  equals  8), 
indicating  the  user  desires  the  Log-Log  transfor 
mation,  ITRAN  is  set  equal  to  3. 

If  PFK  9  has  been  depressed  (NDET(4)  equals  9), 
indicating  the  user  desires  the  Arc  Sin  transfor 
mation,  ITRAN  is  set  equal  to  4. 

If  PFK  10  has  been  depressed  (NDET(4)  equals  10) 
indicating  the  user  desires  the  Weibull  transfor 
mation,  ITRAN  is  set  equal  to  5. 

If  any  other  programmed  function  key  has  been 
depressed,  return  to  detain  statement  2. 

KEY  is  set  equal  to  50 
KEY  is  set  equal  to  1. 

Display  variable  DV1  is  re-initialized. 

The  screen  is  cleared  of  all  existing  displays. 


23 
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SUBROUTINE  CALC 

Purpose:  This  is  a  conventional  maximum-likelihood  routine  for  the 
quantal  analysis. 

Output : 

P  -  a  single-dimensioned  field,  in  floating  point,  which 

contains  the  observed  proportion  of  successes  for  each 
batch 

PCAP  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  predicted  proportion  of  successes  for  each 
batch 

YPR  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  predicted  y-scores  for  each  batch 

Y  -  a  single-dimensioned  field,  in  floating  point  contain¬ 
ing  the  observed  y-scores  for  each  batch 

YWORK  -  a  single- dimensioned  field,  in  floating  point,  contain¬ 
ing  the  "working"  quantities  of  each  batch 

W  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  weights  of  each  batch 

ER  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  error  of  prediction  of  the  y-scores  in  each 
batch 

SSR  -  a  floating  point  number  which  is  the  sum-of -squares 
for  regression 

SSE  -  a  floating  point  number  which  is  the  sum-of-squares  for 
error 

ALF  -  a  floating  point  number  which  estimates  a 

BETA  -  a  floating  point  number  which  estimates  B 

The  output  variables  are  stored  on  a  direct  access  disk,  logical  unit  28. 


Procedure: 


54 


FORTRAN 

STATEMENT  .  STATEMENT 
NUMBER  NUMBER 

5 


6 


7-8 


9  101 


10-11  3 


12-19 


20-32  77 


COMMENT 

KMAX  is  set  equal  to  0,  indicating  termination 
when  alpha  or  beta  converge*. 

K,  the  number  of  the  cycle,  is  initialized 
equal  to  0. 

ALFO  and  BETO,  the  initial  estimates  of  ALF  and 
SETA  are  initialized  equal  to  0. 

K,  the  number  of  the  cycle  is  tested.  If  K 
is  negative,  go  out;  if  K  is  equal  to  zero,  pro¬ 
ceed  to  start;  and  if  K  is  positive,  proceed  to 
FORTRAN  statement  21  (iteration) . 

Loop  58  -  The  number  of  successes  are  translated 
into  proportion  of  successes  in  each  class. 

Loop  5  -  To  avoid  singularities,  if  a  class  has 
no  successes,  it  will  be  translated  as  though 
it  had  H  successes.  If  all  individuals  in  some 
class  show  success,  (n  -  H)  successes  will  be 
recorded. 

Loop  11  -  Y(Observed)  is  calculated  according 
to  the  chosen  transformation:  for  the  Probit 
transformation,  function  YORMP  is  called;  for 
the  Logit  -  DFP;  for  the  Log-Log  -  DFP3;  for 


*This  subroutine  was  adopted  from  a  batch-mode  program  which  permitted 
other  modes  of  termination. 


STATE^NT 

NUMBER 

33 


34-36 

37 

38-52 

53-54 

55 

56 


SUBROUTINE  CALC  (Continued) 

FORTRAN 

STATEMENT  COMMENT 
NUMBER 


55 


the  Arc  Sine  transformation  -  DFP4;  and  for 
the  Weibull  transformation,  DFP2  is  called. 

K,  the  number  of  cycles  is  tested.  If  K  is 
negative,  go  out;  if  K  is  equal  to  zero,  proceed 
to  FORTRAN  statement  10;  and  if  K  is  positive, 
proceed  to  FORTRAN  statement  21. 

10  Loop  12  -  Initially  the  working  y-scores  are 

set  equal  to  the  observed  y-scores  and  the 
weights  are  set  equal  to  the  n^. 

Bypass  recalculation  if  this  is  initial  loop. 

21  Loop  23  -  The  Z-field  is  calculated  according 

to  the  chosen  transformation:  for  the  Probit 
transformation,  function  YORMZ  is  used;  for  the 
Logit  transformation,  DFZ  is  used;  for  the  Log- 
Log  transformation,  DFZ3  is  used;  for  the  Arc- 
Sin  transformation,  DFZ4  is  used;  and  for  the 
Weibull  transformation,  DFZ2  is  used. 

38  No  function  here  (came  from  batch  program) . 

41  If  the  cycle  number  is  negative  (error),  go  out; 

otherwise,  proceed  to  FORTRAN  statement  84. 

84  The  weighted  regression  subroutine,  iv'REG,  is 

called  to  estimate  the  predicted  y-scores,  alpha. 


and  beta. 


SUBROUTINE  CALC  (Continued) 
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STATEMENT 

FORTRAN 

STATEMENT 

COMMENT 

NUMBER 

NUMBER 

57 

. 

The  number  of  the  cycle  is  increased  by  one. 

58. 

If  ALF  is  relatively  close  to  ALFO  or  BETA  is 

60 

60 

relatively  close  to  BETO,  proceed  to  the  next 

statement.  Otherwise,  proceed  to  FORTRAN  state¬ 
ment  60. 

ALFO  is  re- initialized  equal  to  ALF. 

61 

BETO  is  re-initialized  equal  to  BETA. 

62 

If  K,  the  cycle  number  is  less  than  5,  proceed 

to  FORTRAN  statement  141;  otherwise,  proceed  to 

FORTRAN  statement  99. 

63-75 

141 

Loop  65  -  PCAP,  the  predicted  proportion  of 

successes  is  calculated  for  each  batch  according 

to  the  chosen  transformation:  for  the  Probit 

transformation,  YORMX  is  used;  for  the  Logit 

transformation,  DFX  is  used;  for  the  Log-Log, 

DFX3  is  used;  for  the  Arc  Sin,  DRX4  is  used;  and 

76 

for  the  Weibull  transformation,  DFX2  is  used. 

Make  another  iteration. 

78-86 

99 

P,  PCAP,  YPR,  Y,  YWORK,  W,  ER,  SSR,  SSE,  ALF, 

87 

2001 

and  BETA  are  written  onto  direct  access  disk 

Control  is  returned  to  the  calling  program. 

PROBIT  FUNCTIONS:  YORMZ(Y) 

YORMX(Y) 

YORMP(P) 

Definitions: 

Z  =  YORMZ(Y)  Z  »  [! / (2TI)*5]  e-^ 

P  =  YORMX(Y)  P  =  [l/(2ir)^]  £  e dt 


57 


Y  =  YORMP(P) 


obtains  Y  from  P  as  defined  in  YORMX 


Procedure: 

YORMZ: 

YORMX: 


trivial 

If  | Y  |  ^3.0,  YORMX  is  evaluated  by  the  continued  fraction 


expansion  .  ,2 

>l  *  [1/(2.)*]  e  *  {  L.  t- 


3+ 


2  2 

2*1 

5-  7+ 


which  is  continued  until  the  terms  begin  declining  and  the 

-14 

relative  error  is  less  than  10 
If  | Y |  >  3.0,  the  continued  fraction 

2 


1  -  [l/(2*n  e 


*1  B-*  {  i_  !_  i_ 


y+  y+  y+  y+ 


...  } 


is  employed. 

YORMP:  As  a  first  guess,  the  Hastings  form  for  inverting  the  erfc- 
function  is  employed.  This  is  improved  by  a  Newton  itera¬ 
tion,  until  the  relative  error  (P  calculated  from  YORMX  - 
P  input)/ (P  input)  is  less  than  10~^. 


PROBIT  FUNCTIONS:  (Continued) 


Limitations: 

If  P  0,  -DOMEG  is  returned;  if  P  _>  1,  +D0MEG  is  returned, 
where  DOMEG  =  .99999999  x  1038. 


LOGIT  FUNCTIONS:  DFZ(Y) 
DFX(Y) 
DFP(P) 


Definitions: 


Z  =  DFZ(Y) 

N 

II 

% 

>—• 

+ 

to 

X 

w 

P  *  DFX(Y) 

P  *  ey/(l  +  ey) 

Y  =  DFP(P) 

Y  »  log  (P/(l  -  P)) 

LOG  -  LOG  FUNCTIONS: 

DFZ3(Y) 

DFX3(Y) 

DFP3(P) 

Definitions: 

Z  =  DFZ3(Y) 

Z>  -P  log(P) t  where  P 

-Y 

_  -e 

P  =  DFX3(Y) 

P  ■  e 

Y  =  DFP3CP) 

Y  =  -log(-Iog(P)) 

DFZ4(Y) 

DFX4(Y) 

DFP4 (P) 

Z  =  h  cos(Y) 

.  p  =  4  (1  +  sin(Y)) 
Y  =  arsin(2P  -  1) 


ARC  SIN  FUNCTIONS: 

Definitions: 

Z  =  DFZ4(Y) 

P  =  DFX4(Y) 

Y  =  DFP4(P) 


WE I BULL  FUNCTIONS:  DRZ2(Y) 

DFX2(Y) 

DFP2(P) 


60 


Z  «  y^V^' 
p  s  i  _  e"Y 

Y  =  [-r  Icg(l  -  P)]1/r 

SUBROUTINE  WREG: 

This  is  a  conventional  weighted  regression  routine  to  obtain  SSR, 

the  sun- of- squares  for  regression;  SSE,  the  sum- of- squares  for  error; 

2 

VSE,  the  mean  squares  for  error;  ALF;  BETA;  and  R  ,  which  is  SS  Reg/SS 
Total.  The  weighted  regression  routine  is  applied  to  the  "working" 
quantities  while  the  determination  of  error  uses  the  observed  quantities 
as  input. 


Definitions: 
Z  *  DFZ2(Y) 

P  *  DFX2(Y) 

Y  =»  DFP2(P) 
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SUBROUTINE  GRACDF 


Purpose:  The  purpose  of  this  subroutine  is  to  display  a  graph  of  the 
observed  and  predicted  proportion  of  successes  for  each  batch 
(class).  • 


Input: 


P  -  a  single-dimensioned  field,  read  from  direct  access  disk, 
containing  the  observed  proportion  of  successes  for  each 
batch. 

PCAP  -  a  single-dimensioned  field,  read  from  direct  access  disk, 
containing  the  predicted  proportion  of  successes  for  each 
batch 


Procedure: 


FORTRAN 

STATEMENT  STATEMENT 

NUMBER  NUMBER 

5 

6-7 

8 

9 

10-14 


\5-3* 


COMMENT 

Display  variable  DVD  is  established. 

P  and  PCAP  are  read  from  direct  access  disk  28. 
Display  variable  DVD  is  re- initialized. 

The  scope  is  cleared  of  all  existing  displays. 
Loop  10  stores  the  integers  1  to  NCL  in  the 
single-dimensioned  field  IS,  field  X  into 
single-dimensioned  field  XS,  field  C  into 
single-dimensioned  field  PS,  and  field  PCAP 
into  single-dimensioned  field  PCAPS. 

Loop  500  sorts  XS  into  ascending  order,  rearrang¬ 
ing  the  number  of  the  batch  (IS),  the  observed 
proportion  of  successes  (PS)  and  the  ;  redicted 
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SUBROUTINE  GRACDF  (Continued) 

FORTRAN 

STATEMENT  STATEMENT  COMMENT 

NUMBER  NUMBER 

proportion  of  successes  (PCAPS)  in  corresponding 
order. 

33  XMIN,  the  minimum  of  the  x-dosages  is  set  equal 

to  XS(1). 

34  XMAX,  the  maximum  of  the  x-values  is  set  equal 

to  XS(NCL) . 

35-36  XUP  and  XLOW  are  calculated  for  the  upper  and 

lower  coordinates  of  the  horizontal  axis  of 
the  graph. 

37-40  FORMAT  statement  1000,  an  explanation  of  the 

graph,  is  plotted  on  the  screen. 

41  The  coordinates  of  the  lower  left  comer  of 

the  screen  (XLOW,  -.2)  and  the  upper  right  cor¬ 
ner  (XUP,  2.)  are  established  by  the  UCORD  com¬ 
mand. 

42-43  The  vertical  axis  is  positioned  on  the  screen. 

44-45  The  horizontal  axis  is  positioned  on  the  screen. 

46-47  A  message  in  FORMAT  statement  110  is  written  for 

systems  programming  use. 

48-55  Loop  200  constructs  a  step  function  of  the  pre¬ 

dicted  proportion  of  successes  for  each  x-dosage 
arranged  in  ascending  order. 
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SUBROUTINE  GRACDF  (Continued) 

FORTRAN 

STATEMENT  STATEMENT  COMMENT 

NUMBER  NUMBER 

56-62  .  Loop  201  places  asterisks  on  the  diagran  to 

denote  the  observed  proportion  of  successes  for 
each  dosage. 

63-69  Loop  16  constructs  the  markings  on  the  vertical 

scale  and  writes  the  proportion  of  successes  in 
tenths  below  this  axis. 

70-82  Loop  17  constructs  the  markings  on  the  horizon¬ 

tal  axis  and  writes  the  number  of  the  batch  for 
each  x-dosage  below  this  axis. 

83-90  Thu  minimum  and  maximum  dosage- values  are  posi¬ 

tioned  on  the  lower  portion  of  the  screen. 

91  The  size  routine  is  employed  to  determine  whether 

or  not  the  orders  of  display  variable  DVD  will 
fit  into  the  space  that  is  available  in  the  Buffer 
Table. 

92  If  the  orders  of  display  variable  DVD  will  not 

fit  into  the  available  space  in  the  BT,  proceed 
to  FORTRAN  Statement  50. 

93  Display  variable  DVD  is  plotted  on  the  screen. 

94  75  Procedure  is  detained  until  the  user  wishes  to 

continue. 

95  If  an  interrupt  has  been  given  by  the  user  from 

any  device  except  the  programmed  function  ke\ - 


! 
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STATEMENT 

NUMBER 


96 


97 


98 


99 

100 
101 
102 

103-104 


105 

106 


FORTRAN 

STATEMENT 

NUMBER 


750 

67 

SO 


74 


SUBROUTINE  GRACDF  (Continued) 


COHMENT 


board  (K  is  not  equal  to  1),  proceed  to  detain 
statement  75. 

If  programed  function  key  (PFK)  31  has  been 
depressed  (NDET(4)  equals  31),  indicating  the 
user  wishes  to  terminate  the  program,  proceed 
to  FORTRAN  statement  750  (exit) . 

■  If  PFK  30  has  been  depressed  (NDET(4)  equals 
30),  indicating  the  user  wishes  to  restart  the 
program,  proceed  to  FORTRAN  statement  67. 

If  PFK  1  has  been  depressed  NDET(4)  equals  1), 
indicating  the  user  wishes  to  proceed  in  the 
program,  proceed  to  FORTRAN  statement  74. 

Return  to  detain  statement  75. 

KEY  is  set  equal  to  50. 

Proceed  (brunch  around  Restart). 

Restart. 

II,  the  size  cf  cpace  neccessary  to  display  the 
display  variable  DVD  is  written  for  the  program¬ 
mer’s  USj. 

The  screen  is  cleared  of  all  existing  displays. 
Display  variable  DVD  is  re-init  alized. 


SUBROUTINE  OUTPUA 


The  purpose  of  the  OUTPUA  subroutine  is  to  display  detailed 
numerical  results  including  the  size  of  each  batch,  the  ob¬ 
served  proportion  of  successes,  the  predicted  proportion  of 
success  under  the  quantal  model,  the  predicted  y-scores  (Pro¬ 
bits,  etc.),  the  observed  y-scores,  the  "working"  quantities, 
the  weights,  the  observed  x-dosages,  and  the  error  of  pre¬ 
diction  of  the  y-scorss  for  each  class.  A  maximum  of  ten 
classes  is  presented  on  the  screen  at  a  time  —  five  classes 
on  the  upper  portion  of  the  screen  and  five  more  classes  on 
the  bottom  half  of  the  screen.  This  section  has  been  divided 
into  "pages"  so  that  the  user  may  return  to  the  first  page 
of  the  output  at  any  time  by  depressing  programmed  function 
key  4. 

P  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  observed  proportion  o*  successes  of  each  batch 
PCAP  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  predicted  proportion  of  successes  in  each  batch 
YPR  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  predicted  y-scores  for  each  batch 
Y  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  observed  y-scores  for  each  batch 
YWORK  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  "working"  quantities  of  each  batch 
W  -  a  single-dimensioned  field,  in  floating  point,  contain¬ 
ing  the  weights  in  each  batch 
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ER  -  a  single- dimensioned  field,  in  floating  point,  contain¬ 
ing  the  error  of  prediction  of  the  y-scores  in  each  batch 
The  above  input  variables  are  stored  on  direct  access  disk  28. 

A  -  a  single-dimensioned  field,  containing  the  names  of  the 
five  transformations  available  to  the  user-  defined  in 
the  DATA  statement. 

Procedure: 


FORTRAN 

STATEMENT  STATEMENT 

NUMBER  NUMBER 

5 

6-12 

13-14 

IS  102 


16 


17 


18 


19 


20 


COMMENT 

Display  variable  DVA  is  established. 

P,  PCAP,  YPR,  Y,  YWORK,  W,  ER  are  read  from  the 
direct  access  disk  28. 

IXII  and  IXI  are  calculated,  to  be  employed 
later  in  data  field- A. 

NL,  indicating  the  number  of  classes  to  be 
displayed,  is  set  equal  to  NCL. 

JKEY  is  set  equal  to  zero.  JKEY  =  1  indicates 
a  completed  page  of  the  calibration  chart. 

JEY  is  set  equal  to  zero.  JEY  =  1  indicates 
that  this  is  the  last  page  of  the  output. 

I,  indicating  the  page  number,  is  initialized 
equal  to  one. 

K  is  set  equal  to  zero.  K  is  an  even  number 
indicating  the  top  portion  of  the  screen. 

NI,  indicating  the  increment  of  the  classes  to 
be  displayed,  is  set  equal  to  one. 
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STATEMENT 

NUMBER 

21 

22 


23 

24-26 

27-29 

30-59 

61-63 

64 


65 

66 

67 

68 


FORTRAN 

STATEMENT 

NUMBER 

94 


1 


3 

6 

11 


9 

100 


SUBROUTINE  OUTPUA  (Continued) 

COMMENT 

If  NL  is  less  than  10,  JEY  is  set  equal  to  1. 

NL  is  the  number  of  classes  left  to  be  displayed. 
If  NL,  the  number  of  classes  left  to  be  display¬ 
ed,  is  greater  than  5,  proceed  to  FORTRAN  state¬ 
ment  2. 

JKEY  =  1  indicates  a  completed  page  of  output. 

NN  controls  the  amount  displayed  ir  each  section 
(two  sections  per  page). 

Decision  whether  one  or  two  sections  are  needed 
on  the  page. 

Output  via  display  variable  DVA. 

Sequencing  of  pages  to  be  displayed. 

If  JKEY  is  equal  to  1,  indicating  completed 
page  of  calibration  chart,  proceed  to  FORTRAN 
statement  74;  otherwise,  proceed  to  FORTRAN 
statement  94. 

Display  variable  DVA  is  plotted  on  the  screen. 
Procedure  is  detained  to  permit  the  use:  to 
change  pages. 

Display  variable  DVA  is  re-initialized.  ' 

If  an  interrupt  has  been  given  by  the  user  from 
any  device  except  the  programmed  function  key¬ 
board  (J  is  not  equal  to  1),  return  to  detain 


statement  100. 
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STATEMENT 

NUMBER 

69 


70 


71 


72 


74 

76 

77 

78-80 


81 


82-84 


SUBROUTINE  OUTPUA  (Continued) 

FORTRAN 

STATEMENT  COMMENT 
NUMBER 

If  programmed  function  key  (PFK)  31  has  been 
depressed  (NDET(4)  equals  31),  indicating 
the  user  wishes  to  terminate  the  program, 
proceed  to  FORTRAN  statement  750  (exit) . 

If  PFK  30  has  been  depressed  (NDET(4)  equals 
30) ,  indicating  the  user  wishes  to  restart 
.  the  program,  proceed  to  FORTRAN  statement  67. 
If  PFK  1  has  been  depressed  (NDET(4)  equals 
1) ,  indicating  the  user  wishes  to  continue  to 
the  next  page  of  the  output,  proceed  to 
FORTRAN  statement  101. 

If  PFK  4  has  been  depressed  (NDET(4)  equals 
4) ,  indicating  the  user  wishes  to  return  to 
page  1  of  the  output,  proceed  to  FORTRAN 
statement  102. 

101  The  scope  is  cleared  of  all  existing  displays 

4  12  -  the  next  page  -  is  equal  to  I  +  1. 

The  scope  is  cleared  of  all  existing  displays 
The  title  message,  in  FORMAT  statement  5,  is 
plotted  on  the  screen. 

If  JEY  equals  1,  indicating  the  last  page  of 
the  output,  proceed  to  FORTRAN  statement  12. 
The  sub-title  message,  in  FORMAT  statement  60 
is  written  onto  display  variable  DVA. 
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SUBROUTINE  OUTFUA  (Continued) 


STATEMENT 

number 

FORTRAN 

STATEMENT 

NUMBER 

COMMENT 

85 

• 

I*  tIle  Pa8e  number,  is  incremented  by  one. 

87-89 

12 

The  sub-title  message,  in  FORMAT  statement 

• 

65  is  written  onto  display  variable  DVA. 

91 

750 

KEY  is  set  equal  to  50. 

93 

67 

KEY  is  set  equal  to  1. 

94 

74 

The  scope  is  cleared  of  all  existing  displays. 

95 

•  Display  variable  DVA  is  re-initialized. 

•  SUBROUTINE  OUTPUT 
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Purpose:  The  purpose  of  the  OUTPUT  subroutine  is  two-fold: 

(1)  display  an  analysis  of  variance  table  and 

(2)  display  the  graph  of  the  straight  line  Y  «  a  +  Bx 

where  a  and  3  have  been  estimated  in  the  analysis. 

Input : 

ALF  •?  a  floating  point  number  obtained  from  the  analysis, 
which  estimat  es  a 

BETA  -  a  floating  point  number  obtained  from  the  analysis, 
which  estimates  p 

SSR  -  a  floating  point  number  obtained  from  the  analysis 
which  is  the  sum-cf-squares  for  regression 

SSE  -  a  floating  point  number  obtained  from  the  analysis 
which  is  the  sum-of-squares  for  error 

Y  -  a  single-dimensioned  field,  in  floating  point 
containing  the  y-scores  for  each  class 

The  above  input  variables  are  obtained  from  a  direct  access  disk. 

A  -  a  single-dimensioned  field,  containing  the  names  of 
the  five  transformations  available  to  the  user  - 
defined  in  the  DATA  statement 

Procedure: 

FORTRAN 

STATEMENT  STATEMENT  COMMENT 

NUMBER  NUMBER 

5  Display  variable  DVB  is  established. 

6  The  value  of  KEY  is  tested  to  determine  if  the 

user  is  returning  to  the  QUANTAL  program  after 
calling  the  CALCG  program  (KEY  =  25).  If 

the  value  of  KEY  is  25,  even  these  values 


SUBROUTINE  OUTPUT  (Continued) 

FORTRAN 

STATEMENT  STATEMENT  COMMENT 

NUMBER  NUMBER 


71 


which  were  in  COMMON  need  to  be  re-read: 


7-10  ‘  The  values  of  X,  N,  NR,  NCL  and  ITRAN  are 

read  from  direct  access  disk  28. 

12  1000  The  field  Y  is  read  from  direct  access  disk  28. 

13-14  2000-  Variables  ALF,  BETA,  SSR,  and  SSE  are  read 

2001 

into  core. 


15-16  IXI  and  IXII  are  calculated  to  obtain  the 

correct  subscript  for  field  A. 

17  JR,  the  degrees  of  freedom  for  regression,  is 

set  equal  to  1. 

18  JE,  the  degrees  of  freedom  for  error,  is  set 

equal  to  NCL- 2. 

19  JTOT,  the  total  degrees  of  freedom  is  calcu¬ 
lated. 

20  On  return  from  dosage-transformation  routine, 

the  analysis  of  variance  display  is  skipped. 

21  VSE,  the  mean  squares  for  error  is  calculated. 

22  ST0T,  the  sum-of-squares  total  is  ca]culated. 

23  The  F-statistic  is  calculated. 

24-30  FORMAT  statements  71  and  72,  displaying  the 

Analysis  of  Variance  table,  are  plotted  on 
the  screen. 


31  73  Procedure  is  detained  until  user  has  read  the 


analysis  of  Variance  table. 
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STATEMENT 

NUMBER 

32 


33 


34 


'  35 


36 


37 

38 

39-41 


SUBROUTINE  OUTPUT  (Continued) 

FORTRAN 

STATEMENT  COMMENT 
NUMBER 

If  an  interrupt  has  been  given  bv  the  user 
from  any  other  device  than  the  programed 
function  keyboard  (I  is  not  equal  to  1) , 
return  to  detain  statement  73. 

If  programmed  function  key  (PFK)  31  has  been 
depressed  by  the  user  (NDET(4)  equals  31), 

•  indicating  the  user  wishes  to  terminate  the 
program,  proceed  to  FORTRAN  statement  750 
(exit) . 

If  PFK  30  has  been  depressed  (NDET(4)  equals 
30) ,  indicating  the  user  wishes  to  restart  the 
program,  proceed  to  FORTRAN  67. 

If  PFK  11  has  been  depressed  (NDET(4)  equals 
11),  indicating  the  user  wishes  to  leave  the 
QUANTAL  program,  proceed  to  FORTRAN  statement 
80. 

If  any  PFK  other  than  key  1  or  those  previous¬ 
ly'  mentioned,  has  been  depressed  (NDET(4)  is 
not  equal  to  1),  return  to  detain  statement  73. 
491  The  scope  is  cleared  of  all  existing  displays. 

Display  variable  DVB  is  re-initialized. 

FORMAT  statement  400,  a  title  message,  is 
written  onto  display  variable  DVB. 
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SUBROUTINE  OUTPUT  (Continued) 


FORTRAN 

STATEMENT 

STATEMENT 

COMMENT 

NUMBER 

NUMBER 

42-44 

•  2002 

The  message  in  FORMAT  statement  2004  is  writ¬ 

ten  onto  display  variable  DVB. 

45-52 

• 

Loop  4  -  The  maximum  of  the  x-values  and  the 

maximum  of  the  y-values  are  determined. 

53-60 

In  loop  8,  the  minimum  of  the  x-values  and  the 

minimum  of  the  y-values  are  determined. 

61-71 

•  XLGTH,  XINC,  YLGTH,  YINC,  Y2MIN,  Y2MAX  and 

X2MAX  are  calculated  to  determine  the  coordi¬ 

nates  of  the  screen. 

72 

The  coordinates  of  the  lower  left  corner  of 

the  screen  (XMIN,  Y2MIN)  and  the  upper  right 

comer  (X2MAX,  Y2MAX)  are  established  by  the 

• 

UCORD  statement. 

73-74 

The  vertical  axis  is  positioned  on  the  screen. 

75-76 

The  horizontal  axis  is  positioned  on  the  screen 

77 

13 

Display  variable  DVB  is  plotted  on  the  screen. 

78-79 

Diagnostic. 

80 

14 

Tie  length  of  the  x-axis  is  divided  into  200 

-  equal  lengths  (XLGTH) . 

81 

If  KEY  equals  6,  indicating  the  x-dosages  have 

been  altered  by  one  of  the  functions,  proceed 

to  FORTRAN  statement  270  and  skip  loop  15. 

74 


STATEMENT 

NUMBER 

82-88 

89-92 

93-94 

97-98 

99-136 


137 

138 
140 


SUBROUTINE  OUTPUT  (Continued) 

FORTRAN 

STATEMENT  COMMENT 

NUMBER 

In  loop  15,  the  points  on  the  straight  line 
Y  =  a  +  8x  are  calculated  and  positioned  onto 
display  variable  DVB. 

270  Loop  16  positions  an  asterisk  onto  display 

variable  DVB  at  each  observed  y-value. 

XX  and  XNUM  are  calculated  to  determine  ten 
•  equally  spaced  x- values  on  the  horizontal  axis. 
36  XX  and  XNUM  are  calculated  to  determine  ten 

equally  spaced  y-values  on  the  vertical  axis. 

26  Loop  24  positions  the  markings  on  the  horizon¬ 

tal  scale  when  JEY  equals  1  and  the  markings 
on  the  vertical  scale  when  JEY  equals  2  at 
ten  equally  spaced  values.  The  x-  and  y- 
values  are  aliened  by  these  markings  by  the 
five  different  FORMAT  statements. 

If  JEY  is  not  equal  to  1,  indicating  both  axes 
have  been  marked,  proceed  to  FORTRAN  statement 
22. 

JEY  is  set  equal  to  2,  to  initiate  y-axis. 

22  The  size  routine  is  employed  to  determine 

whether  or  not  the  orders  of  display  variable 
DVB  will  fit  into  the  space  that  is  available 


in  the  BT. 


75 


SUBROUTINE  OUTPUT  (Continued) 


FORTRAN 

STATEMENT 

STATEMENT 

COMMENT 

NUMBER 

NUMBER 

141 

If  the  orders  of  display  variable  DVB  will 

* 

not  fit  into  the  available  space  in  the  BT, 

• 

proceed  to  FORTRAN  statement  50. 

142 

Display  variable  DVB  is  plotted  on  the  screen. 

143-144 

50 

Diagnostic. 

14S 

75 

Procedure  is  detained  until  the  user  has  read 

•  the  display. 

146 

If  an  interrupt  has  been  given  by  the  user 

from  any  other  device  except  the  programmed 

function  keyboard  (K  is  not  equal  to  1), 

return  to  detain  .statement  75. 

147 

If  PFK  31  has  been  depressed  (NDET(4)  equals 

• 

31),  indicating  the  user  wishes  to  terminate 

the  program,  proceed  to  FORTRAN  statement  750 

(exit). 

148 

If  PFK  30  has  been  depressed  (NDET(4)  equals 

30),  indicating  the  user  wishes  to  restart  the 

program,  proceed  to  FORTRAN  statement  67. 

149 

If  PFK  1  has  been  depressed  (NDET(4)  equals 

1) ,  indicating  the  user  wishes  to  proceed  in 

the  program,  exit  from  the  subroutine. 

151-154 

80 

X,  N,  NR,  NCL  and  ITRAN  are  written  onto 

direct  access  disk  28. 


SUBROUTINE  OUTPUT  (Continued) 
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FORTRAN 

STATEMENT  STATEMENT  COMMENT 


NUMBER 

NUMBER  - 

1SS 

• 

KEY  =  50  to  leave  QUANTAL  in  order  to  do 

calculations. 

1S7 

750 

Prepare  for  termination. 

159 

67 

KEY  is  set  equal  to  1  (user  wants  to  restart) . 

160 

99 

The  screen  is  blanked  of  all  existing  displays. 

161 

Display  variable  DVB  is  re-initialized. 
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SUBROUTINE  CHNGFN 

Purpose:  The  purpose  of  the  CHNGFN  subroutine  is  to  make  various  func¬ 
tions  available  for  possible  transformations  of  the  x-dosages. 

Procedure: 

FORTRAN 

STATEMENT  STATEMENT  COMMENT 
NUMBER  NUMBER 

3  The  display  variables  DVA,  DVB,  and  DVE  are 

established. 

4  Display  variable  DVB  is  re-initialized. 

5  The  scope  is  cleared  of  all  existing  displays. 

6-9  FORMAT  statement  5  is  plotted  on  the  screen. 

10  1  Procedure  is  detained  until  the  user  decides 

whether  or  not  he  wishes  to  change  the  x-dos- 
ages . 

11  If  an  interrupt  has  been  given  by  the  user 

from  any  device  except  the  programmed  function 
keyboard  (J  is  not  equal  to  1),  return  to 
detain  statement  1. 

12  If  programmed  function  key  (PFK)  31  has  been 

depressed  (NDET(4)  equals  31),  indicating 

the  user  wishes  to  terminate  the  program,  pro¬ 
ceed  to  FORTRAN  statement  750. 

13  If  PFK  30  has  been  depressed  (NDET(4)  equals  30), 

indicating  the  user  wishes  to  restart  the 
program,  proceed  to  FORTRAN  statement  99 

(exit) . 


SUBROUTINE  CHNGFN  (Continued) 
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STATEMENT 

NUMBER 

14 


15 


17-20 


21 


22 


23 


24 


FORTRAN 

STATEMENT  COMMENT 
NUMBER 


If  PFK  1  has  been  depressed  (NDET(4)  equals 
1),  indicating  the  user  wishes  no  alteration 
to  be  made  on  the  x- dosages,  proceed  to 
FORTRAN  statement  50. 

If  PFK  3  has  been  depressed  (NDET(4)  equals 
3),  indicating  the  user  desires  to  alter  the 
x-dosages  by  one  of  the  functions,  ‘proceed  to 
FORTRAN  statement  10. 

10  The  statement  to  obtain  which  function  the 

user  wishes  (FORMAT  2)  is  plotted  on  the  screen. 
A  cursor  is  placed  in  the  first  unprotected 
position  of  FORMAT  statement  2  so  that  the 
user  can  type  in  the  function  desired. 

3  Procedure  is  detained  until  the  user  types  in 

the  appropriate  function. 

If  PFK  31  has  been  depressed,  (NDET(4)  equals 
31),  indicating  the  user  wishes  to  terminate 
the  program,  proceed  to  FORTRAN  statement 
750  (exit). 

If  PFK  30  has  been  depressed  (NDET(4)  equals 
30),  indicating  the  user  wishes  to  restart 
the  program,  proceed  to  FORTRAN  statement  99. 


SUBROUTINE  CHNGFN . (Continued) 
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STATEMENT 

NUMBER 

25 


26-28 

29 

30 

31 


32 


'  33 

35 

36 

37-40 


FORTRAN 

STATEMENT-  COMMENT 
NUMBER 

If  an  interrupt  has  been  given  by  the  user 
from  any  device  except  the  alphameric  type¬ 
writer  (J  is  not  equal  to  2) ,  return  to  de¬ 
tain  statement  ?. 

XFUN,  the  real  value  for  IFUN  is  read. 

,  The  pointers  for  dummy  unit  4  are  reset. 

IFUN  =  XFUN. 

If  the  number  of  the  function  desired  to 
alter  the  x-dosages  is  less  than  zero  or 
greater  than  12,  proceed  to  FORTRAN  state¬ 
ment  20  to  write  an  error  message. 

When  IFUN  has  the  value  8,  9,  or  10  proceed 
to  FORTRAN  statement  1050  to  ascertain  that 
all  x-dosages  are  non-negative  real-numbers. 
If  IFUN  has  any  other  value  (less  than  12), 
proceed  to  FORTRAN  statement  53. 

53  KEY  is  set  equal  to  6,  indicating  subroutine 

FNS  should  be  called. 

20  Incorrectly  answered  FORMAT  statement  2  is 

cleared  from  the  screen. 

Display  variable  DVA  is  re-initialized. 

The  error  message  in  FORMAT  6  is  plotted  on 


the  screen. 


SUBROUTINE  CHNGFN  (Continued) 
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FORTRAN 

STATEMENT 

STATEMENT 

COMMENT 

NUMBER 

NUMBER 

41 

• 

Display  variable  DVE  is  re-initialized. 

43-45 

1050 

Loop  55  tests  all  x-values  to  ascertain  they 

are  non-negative  real  numbers.  If  any  x-value 

is  negative,  proceed  to  write  an  error  message. 

47 

52 

Display  variable  DVA  is  re-initialized. 

48 

The  current  instance  of  DVA  is  removed  from 

the  Buffer  Table. 

49-52 

The  message  in  FORMAT  statement  51  which  states 

that  the  chosen  function  cannot  be  employed 

is  plotted  on  the  screen. 

53 

Display  variable  DVA  is  re- initialized. 

55 

750 

KEY  is  set  equal  to  50. 

57 

99 

KEY  is  set  equal  to  1. 

59 

50 

KEY  is  set  equal  to  5. 

60 

100 

The  scope  is  cleared  of  all  existing  displays. 

61 

Display  variables  DVA,  DVB,  and  DVE  are  re-ini¬ 
tialized. 

Purpose: 


Procedure : 


81 


SUBROUTINE  FNS  . 

The  purpose  of  subroutine  FNS  is  to  transform  the  x-dosages 
by  one  of  twelve  mathematical  functions  made  available  to 
the  user  in  CHNGFN. 

Loop  4  is  executed  for  each  x-dosage.  The  function  employed 


is  determined  by  the  value  of  IFUN  (obtained  in  subroutine 
CHNGFN)  in  the  computed  GO  TO  statement. 
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FIGURE  4.1 

fROGRAM  CONTROL  CARDS 

//QUANTAL  JOB  (.'ACCOUNT#’  ,10,4) ,  'PROGRAMMER.'  ,MSGLEVEL=1 
//STEP1  EXEC  FORTGC 

//FORT.SYSLIN  DD  DSNAME= §§CHAIN (ROOT) ,SPACE=(TRK, (150,10,5)) , 
//  UN IT=SYSDA, DISP= (NEW, PASS) 

//FORT.SYSIN  DD  * 

(source  deck  for  MAIN) 

/* 


//STEP 2  EXEC  FORTGC 

//FORT.SYSLIN  DD  DSNAME=5§CHAIN(LINKA) ,DISP= (MOD, PASS) ,UNIT=SYSDA 
//FORT.SYSIN  DD  * 

(source  deck  for  DATA) 

/* 


//STEP 3  EXEC  FORTGC 

’  //FORT.SYSLIN  DD  DSNAME=  fijCHAIN (LINKD) , DI SP=  (MOD , PASS) , UN IT=SYSDA 
//FORT.SYSIN  DD  * 

(source  deck  for  INPUTA) 

/* 


FIGURE  4.1  (Continued) 


//STEP4  EXEC  FORTGC  _ 

//FORT.SYSLIN  DD  DSNAME=$5CHAIN(LINKE) ,DISP= (MOD, PASS) ,UNIT=SYSDA 
//FORT. SYS IN  DD  * 

(source  deck  for  INPUTB) 

/* 


//STEP5  EXEC  FORrGC 

//FORT.SYSLIN  DD  DSNAME=§§CHAIN(LINKI) ,DISP= (MOD, PASS) ,UNIT=SYSDA 
//FORT. SYS IN  DD  * 

(source  deck  for  CHNGTR) 

/* 


//STEP6  EXEC  FORTGC 

//FORT.SYSLIN  DD  DSNAME=§§CHAIN(LINKB) ,DISP= (MOD, PASS), UNIT=SYSDA 
//FORT. SYSIN  DD  * 

(source  deck  for  CALC) 

b 


//STEP7  EXEC  FORTGC 

//FORT.SYSLIN  DD  DSNAME=5$CHAIN(LINKJ) ,DISP= (MOD, PASS) ,UNIT=SYSDA 
//FORT. SYS IN  DD  * 

(source  deck  for  GRACDF) 


FIGURE  4.1  (Continued) 


//STEP8  EXEC  FORTGC 

//FORT.SYSLIN  DD  DSNAME=§$CHAIN(LINKF) ,DISP= (MOD, PASS) ,UNIT=SYSDA 
//FORT.SYSIN  DD  * 

(source  deck  for  OUTPUA) 

/* 


//STEP9  EXEC  FORTGC 

//FORT.SYSLIN  DD  DSNAME=§SCHAIN(LINKC) ,DISP= (MOD, PASS) ,UNIT=SYSDA 
//FORT.SYSIN  DD  * 

(source  deck  for  OUTPUT) 

/* 


//STEP 10  EXEC  FORTGC 

//FORT.SYSLIN  DD  DSNAME={j§CHAIN (LINKG) , DISP= (MOD,PASS) ,UNIT=SYSDA 
//FORT.SYSIN  DD  * 

(source  deck  for  CHNGFN) 

/* 


//STEP 11  EXEC  FORTGC 

//FORT.SYSLIN  DD  DSNAME=$SCHAIN(LINKH),DISP=(MOD,PASS),UNIT=SYSDA 
//FORT.SYSIN  DD  * 

(source  deck  for  Fr«'S) 
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FIGURE  4.1  (Continued) 

//STEP 12  EXEC  LKED,PARM= (LET, LIST,OVLY, XREF) 

//LKED.SYSLMOD  DD  DSN=SYS1 . GRAPHLIB (QUANTAL) ,DISP=SHR, 

//  SPACE* (TRK ,(0,0)) 

(GRAPHLIB (QUANTAL)  is  to  receive  the  load  module) 

//LKED.SYSLIB  DD  DSN=SYS1.GRAFLIB,DISP=SHR 
(GRAFLIB  contains  graphics  routines) 

//DD  DSN=SYS1. FORTLIB,DISP=SHR 
//DD  DSN=SYS1. GRAPHLIB, DISP=SHR 
//DD  DSN*SYS1 . UGALIB , DISP=SHR 
//DD  DSN=SYS1. LINKLIB,DISP=SHR 

(utility  routines  at  the  University  of  Georgia  system) 

//LKED. MODULE  DD  DSN=$5CHAIN,DISP=0LD 

(55CHAIN  contains  the  object  modules,  as  prepared  in  the  FORT  step; 
they  are  referred  to  under  the  DD  name  MODULE) 

//LKED.SYSIN  DD  * 

INCLUDE  MODULE (ROOT)  (contains  MAIN  program) 

INCLUDE  SYSLIB(IBCOM#) 

INCLUDE  SYSLIB(ARITH#) 

INCLUDE  SYSLIB(FIOCS#) 

INCLUDE  SYSLIB(ADCON#) 

INCLUDE  SYSLIB(IHCUATBL) 

INCLUDE  SYSLIB(IHCUOPT) 

INCLUDE  SYSLIB(ERRMON) 

INCLUDE  SYSLIB(IHCTRCH) 

INCLUDE  SYSLIB(GAFERR) 


FIGURE  4.1  (Continued) 


OVERLAY  ONE 
INCLUDE  SYSLIB(ARCOS) 

INCLUDE  SYSLIB(ATAN2) 

INCLUDE  SYSLIB(COS) 

INCLUDE  SYSLIB(DCOS) 

INCLUDE  SYSLIB (COTAN) 

INCLUDE  SYSLIB (DARCOS) 

INCLUDE  SYSLIB (DEXP) 

INCLUDE  SYSLIB (DLOGIO) 

INCLUDE  SYSLIB(DSQRT) 

INCLUDE  SYSLIB (SQRT) 

INCLUDE  SYSLIB (ARSIN) 

INCLUDE  SYSLIB (ATAN) 

INCLUDE  SYSLIB(SIN) 

INCLUDE  SYSLIB (DSIN) 

INCLUDE  SYS LIB (TAN) 

INCLUDE  SYSLIB(DARSIN) 

INCLUDE  SYSLIB (DLOG) 

OVERLAY  TWOA 

INCLUDE  MODULE (LINKB)  (contains  SUBROUTINE  CALC) 

OVERLAY  TWOA 

INCLUDE  MODULE (LINKH)  (contains  SUBROUTINE  FNS) 

OVERLAY  ONE 

INCLUDE  SYSLIB (DISPLA) 

INCLUDE  SYSLIB (PLOT) 


FIGURE  4.1  (Continued) 
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INCLUDE  SYSLIB (DETAIN) 

INCLUDE  SYSLIB (BUFRS) 

INCLUDE  SYSLIB (CUR$$) 

INCLUDE  SYSLIB (LIGHTS) 

INCLUDE  SYSLIB (SCTDV) 

INCLUDE  SYSLIB (ERASE) 

INCLUDE  SYSLIB (RCUR$) 

INCLUDE  SYSLIB (RESET) 

INCLUDE  SYSLIB($$$$BT) 

INCLUDE  SYSLIB (CHAR) 

INCLUDE  SYSLIB (DETEKT) 

INCLUDE  SYSLIB (READSC) 

INCLUDE  SYSLIB (CLOSE) 

INCLUDE  SYSLIB ($VOVER) 

INCLUDE  SYSLIB ($$OVER) 

INCLUDE  SYSLIB (DUMMY$) 

INCLUDE  SYSLIB(UNPLOT) 

INCLUDE  SYSLIB (PLACES) 

INCLUDE  SYSLIB (BLANK) 

INCLUDE  SYSLIB($ADD$) 

INCLUDE  SYSLIB($$INIT) 

INCLUDE  SYSLIB (SCNDVDK) 

INCLUDE  SYSLIB (WRFMTJ) 

OVERLAY  TWO 

INCLUDE  MODULE(LINKA)  (contains  SUBROUTINE  DATA) 
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OVERLAY  TWO 
INCLUDE  MODULE (LINKD) 
OVERLAY  TWO 
INCLUDE  MODULE (LINKE) 
OVERLAY  TWO 
INCLUDE  MODULE (LINKF) 
OVERLAY  TWO 
INCLUDE  MODULE (LINKG) 
OVERLAY  TWO 

INCLUDE  SYSLIB(CORDCALL) 
INCLUDE  SYSLIB(LINE$$) 
INCLUDE  SYSLIB(SIZE) 
INCLUDE  SYSLIB($CORD$) 
INCLUDE  SYSLIB(POINT$) 
OVERLAY  THREE 
INCLUDE  MODULE (LINKC) 
OVERLAY  THREE 
INCLUDE  MODULE (LINKI) 
OVERLAY  THREE 
INCLUDE  MODULE (LINKJ) 

// 


FIGURE  4.1  (Continued) 


(contains  SUBROUTINE  INPUTA) 
(contains  SUBROUTINE  INPUTB) 
(contains  SUBROUTINE  OUTPUA) 
(contains  SUBROUTINE  CHNGFN) 


(contains  SUBROUTINE  OUTPUT) 

(contains  SUBROUTINE  CHNGTR) 

(contains  SUBROUTINE  GRACDF) 


GRAF  GRAPHIC  ROUTINES 
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ccoi  OOUOLE  PRECISION  FUNCTION  YORMXIOZI 

CCC2  IMPLICIT  RE*L*S(OI 

CC3J  OOUBLE  PRECISION  YORHX 

CCS*  DPI  «  .3989422E040I433 
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ccal  ROUBLE  PRECISION  FUNCTION  0FP4ICPI 

“02  COUBLE  PRECISION  0FP4.DP 

CCCJ  DEP4  *  0PRSIM2.00  •  DP  -  1.00) 

CCC4  RETURN 

CCOS  ENC 
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DOUBLE  PRECISION  FUNCTION  DFP2(DP) 

DOUBLE  PRECISION  DFP2,  DP,  DOMEG,  DRR,  DEX,  DLOG 
COMMON  NDET(5) ,  X(100),  N(100),  NR(IOO),  DV1,  DV2,  DV3, 
1DV4 ,  DV5,  DV6,  0V7,  DVD,  DVE,  DVA,  DVB,  NCL,  KEY, 

2KK,  I FUN,  ITRAN,  IREC 
DOMEG  =  0.9999999999D58 
DFP2  =  DOMEG 
IF(DP.GE.l.DO)  GO  TO  99 
DFP2  =  -DOMEG 
IF(DP.LE.O.DO)  GO  TO  99 
DRR  =  NR(IOO) 

DEX  =  l.DO/DRR 

DFP2  =  (-DRR*DLOG ( 1 . DO-DP) ) **DEX 
99  RETURN 
END 
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0008 


DOUBLE  PRECISION  FUNCTION  DFZ2 (DY) 

DOUBLE  PRECISION  DFZ2 ,  DY,  DR,  DEXP 

COMMON  NDET(5) ,  X(100),  N(100),  NR(IOO),  DV1,  DV2,  DV3, 
1DV4,  DV5,  DV6,  DV7,  DVD,  DVE,  DVA,  DVB,  NCL,  KEY, 

2KK,  I FUN,  ITRAN,  IREC 
LR  =  NR(IOO) 

DR  a  LR 

DFZ2  =  DY** (LR-1)*DEXP(-DY**LR/DR) 

RETURN 

END 


0001  DOUBLE  PRECISION  FUNCTION  DFX2(DY) 

0002  DOUBLE  PRECISION  DFX2,  DY,  DRR,  DEXP 

0003  COMMON  NDET(5) ,  X(100),  N(100),  NR(IOO),  DV1,  DV2, 

1DV3,  DV4,  DV5,  DV6,  DV7,  DVD,  DVE,  DVA,  DVB,  NCL, 

2 KEY,  KK,  I FUN,  ITRAN,  IREC 
0004  DRR  =  NR(IOO) 

0005  DFX2  =  l.DO  -  DEXP (-DY**NR(100)/DRR) 

0006  RETURN 

0007  END 
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CC4S  CALL  L!NE|OVC»>WIN.U) 

CC46  Wft IT£ (6*  110)  IBV.F 

CCW  1X0  FQPHATtlOX,*  PL'FFFR  SPACE  REMAINCNG 

CCA8  CALL  PLACE  t  UVD  »XMlN*0« ) 
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FORTRAN  IV  C  LEV^L  \% 


OUTPUT  DATE  *  7118P  10/48/71 


0001 

0002 

0003 

0004 

0005 

0004 

OOC7 

0008 

OOOS 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0010 

0070 

0021 

0022 

0073 

0024 

007* 

0*26 


0027 

0028 

0079 


0030 

0031 

0037 

0017 

0014 

0075 

0036 

0077 

0038 

0039 

00*0 
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SUBROUTINE  OUTPUT 
01MFNSICN  Allot*  VI 1001 

COMMON  NDCTISI.XUOOt.NI 1001, NR (100 1 ,DVl,DV2,DV3,DV4  .0V5.0V6 .OVT , 
1CVC,0VF,CVA,CVP, NCL,K£Y,KK,IFyN,  IT  RAN, IREC 
OATA  STAR/lH*/,0/lHl/,S/lH-/,A/4H  PR.4H0B1T.4H  l,4H00IT,4H  LOG, 
14H-L 0G.4H  ARC, 4H  S|N,4H  ME  1 , 4HBULL/ 

,  CALL  DISPLAIOVPI 

IF(KEY.NE.75tGO  TO  1000 

RE  AO  (28*  10IX 

Rf  AO  I  28.'  11  IN 

RE A0I2R* 12INR 

RE  AO  1 28'  17  t  NCL  ,  l  TRAN 

KEY  •  0 

1000  RCAOI79'4IV 

2000  °EA0(2R,4|ALFIBETA 

2001  REA0<?P'8)SSR,SSE 

I X 1  «  7  *  ITRAN  -  1 
IX 1 1  «  IXt  ♦  1 
JR  «  1 
JE  «  NCL  -2 
JTOT  »  JR  ♦  JE 
(F(KFY.EQ.6)fn  TO  491 
VSE  *  SSE  /  FLCATIJEI 
STOT  «  SSR  ♦  SSF 
F  •  SSR  /  VSE 

MR  IT  El  4, 71  |A(  1X1  I .  A  (  1X11  I,  SSR, JR, SSR, SSF, JE, VSE,  STOT, JTOT, F 
CALL  WRFMTiirvPI 

71  Fn»M#TCPt'/'Pf*,27X,7A4,»  ANAL VSI S* / • PO*  , 73 X . ' ANAL  VS  I S  OF  VAPI4NC 
5  F  TABLE* /'pn*  ,  21  X  , '  S  UM  OF  * ,  1 1  X,  '  OEGB  EES  OF*  ,  1 2  X,  •  ME  A  N*  /  •  P  SflURCF 
6*,12X.  'SQUARES*  ,12X,*FRFFCCM*  ,  1 3X  ,*  SQUARE  •/•  Pf  REGRESS  t  ON*  ,5X,EiR.« 
7,10X,I5,5X,E18.8/*P0FRRnP* ,12X,et6.B,10<,  1 5 , 7X , F 16. 8 / • PG TOTAL ■ .10X 
R, F18. 8, 10X, 15/* P O', 21X, ‘CALCULATED  F  »  '.E16.8I 

WRITEI4.72IJF 
CAU  WRrMTK  fVBt 

72  FOPMATI  *P0  '/ *P(|  IF  YOU  V I SM  TO  LEAVE  THE  CUANTAl  PROGRAM  TO  OE 

l TERM  INF  F 1 1 , • , I? , ' t , BY*/ *P  CALLING  THE  "CALCG"  PROGRAM,  DEPRESS  PR 
20GPAM  FUNCTION  KEY  11.  IF  YOU'/’P  WISH  TO  CONTINUE,  OEPRFSS  PROG* 
3AM  FUNCTION  KEY  1.* ) 

IBUF  «  PLOTIfVP! 

77  1  «  CETAIN(NCET) 

IFlI.NF.HGn  fC  73 
IFIN0ETI4l.EC.3l  ICC  TC  T5P 
1F(NOET(4|.EC.30)GO  TO  67 
IF (NPET(4».EC.ll IGC  TC  RC 
IF1N0ETI4I.NF.UG0  TO  73 
491  CALL  BLANK 

CALL  RESETIOVBt 

MR  IT  El  4, 430  I Al IX | I, A  I  1X1  I  1 .AIF.RETA 
CALL  MRFMTAICVRI 

40C  FORMAT  1 *P1 •/ *P0  THE  FCLLOMING  IS  THE  GRAPH  OF  TH1  STRAIGHT  LIN 

IE  Y  •  ALPHA  4  BETA  «  XV*P  WHERE  ALPHA  ANO  BETA  HAVE  BEEN  ESTIMATE 
20  IN  THF  '  ,  2  A4 ,  •  ANALYSIS. */*P  HERE  THE  EQUATION  IS  Y  •  *,fl6.6, 
3*  4  ' , F I  6. 6 , '  •  X.*/*P  THE  HO* I ZONT  Al  AXIS  IS  THE  OOSAGE-AXIS  ANO 
4THE  VERTICAL  AXIS  IS  THE  r-AXIS'/'P  the  VALUES  denoteo  BY  •••••  AR 
5f  THE  OBSEPVEn  Y-VALUES-1 

20C2  MR  I T  E ( 4 , 20P4 | 


0042 


123 


FORTRAN  IV  C  LEVEL  !9 


OUTPUT 


0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 

MU 

0112 

Oil? 

0114 

0U3 

0116 

0117 

one 

0119 
0120 
0121 
0122 
0123 
0124 
C  1 25 
0126 
0127 
0128 
0129 
0110 
0131 
0132 
013? 
0134 
0135 
0136 
0137 
0138 
'  0139 
0140 
0141 
0142 
C143 
0144 
0145 

one 

0147 

0148 

0149 

0150 

0151 

0132 

0153 


•  n  *  ,YMiX  "  *  . 

26  cn  ?4  1=1,10 

XX  «  XX  ♦  XNUN 
CO  TO  1 27,28  I , jfy 

27  CALL  PLACFlOVB.XX.o.S 
{ALL  CHARlOVf.O, 1,691 
CALL  PLACE1DVB.XX, VINCI 

GO  TO  29 

28  CALL  PLACE(OVR,0.,XXI 
CALL  CHARI CVt , S, 1 ,69 * 

PlAC£'nVB,*iNC,XX| 

29  f!>f>'*tE..1000.|G0  TO  30 
IF  (XX.  L F. - 100,  1G0  TO  31 
IF (XX. IE .-10. ICO  TO  37 
IFlXX.LT.n. ICO  TO  33 
IFIXX.LT.  10.  IGfl  TO  34 
IFlXX.LT.lOo. ICO  TO  33 
IFIXX.LT.  10^0.  ICC  TO  32 
IF1XX.LT. 10000. ICO  TO  31 

30  WRITFC4.35IXX 
call  wrf h ti i ove i 

35  F00H4T('T.«,Fq.,, 

GO  TO  24 

31  WRITE!*, 40IXX 
CALL  WRFNTilCVPI 

40  F0RNAT(«P..,F7.?| 

GO  TO  24 

32  WRITEC4.451XX 
CALL  WRFHTIIOVBI 

*5  FOPMATI'P.  • ,  F6.2I 
GO  TC  24 

33  WRITEI4,55IX> 
call  wrfmt  1 1 eve  I 

35  H'R^AT I*  P.  .  ,F5.2  I 
GO  TO  24 

34  NR  I TE ( 4 , 60 1 X  X 
r»'.L  WRF«T*(CVR| 

60  F0«NAT(«p.t,F4.2| 

24  CONTINUF 

IF (JEV.NE. 1 IGO  TO  22 
JE  V  .2 
GO  TO  36 

22  !!  *  SI2E1CV6I 

*F 1 1 t.GT. I8UFIG0  TO  50 
IBUF  »  PLOT ( Cv R | 

50  NR  I TE ( 6 , 300 | II 
300  FORMAT  1 10X ,'11  .  .ll5, 

T5  X  •  OFTA ININCET I 
IF1K.NE.UGP  TO  75 
JFIN0ETI4I.E0.31  IGO  TO  ~50 
IF1NOETI4I.EO.30  IGO  TO  67 
IFINDE T14l.EQ.il GO  TO  99 
GO  TO  75 

80  VR  ITEMS' 10IX 
WRITEMR'Uin 
NR ITEI 28' 12 1  nr 


DATE  •  71188 


10/4S/21 
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FORTRAN 

tv  r.  LEVI 

1°  OUTPUT 

015A 

• 

l»P|TF(?«'imCL,  tTPAN 

r.js* 

KFV  x  50 

9156 

CO  TO  99 

01*57 

750 

KPV  «  5r 

rus 

GO  TO  oo 

rjs<j 

.57 

Krv  <  t 

0 1  AO 

99 

CALI  BLANK 

01M 

CAIL  RPSMIOVBI 

"If? 

BP  TURN  . 

01  A? 

ENC 

lo/4R/i  1 


cc«.i  call  RfSfrfovei 


